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bearing  pitch  diameter  (mm) 
ball  diameter  (mm) 

. p 

elastic  modulus  (N/m  ) 

p 

convection  factor  = pc  V z /2  k Ay 
complete  elliptic  integral  of  the  first  kind 
complete  elliptic  integral  of  the  second  kind 
external  force  on  a bearing  (N) 
curvature  factors 

f^-D|j  is  the  radius  of  inner  race  curvature 
f g • Djj  is  the  radius  of  outer  race  curvature 
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m mass  (kg) 

p pressure  (Pa) 
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viscosity-pressure  factor  (1/Pa) 

P viscosity-temperature  factor  (1/°C) 

5 deflection  (mm) 

<Spj^  axial  preload  (mm) 

5p^  radial  preload  (mm) 

T]  lubricant  dynamic  viscosity  (N-sec/m^) 

li  ellipticity  factor,  k = a/b 

/I  coefficient  of  friction  (or  called  traction 

coefficient  for  elastohydrodynamic  lubrication) 

T shear  stress  in  elastohydrodynamic  film  (N/m^) 
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inner  race  rotating  speed  (rad/sec) 
ball  orbiting  speed  about  the  bearing  center 
(rad/sec) 

Poisson's  ratio 

ball  rotating  speeds  about  its  own  axis  (rad/sec) 
spinning  component  of  the  relative  ball/race 
rotating  velocity  (rad/sec) 
rolling  component  of  the  relative  ball/race 
rotating  velocity  (rad/sec) 
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to 
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to 
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to 
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to 
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to 
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ball 

ball/ inner  race  contact 
inner  race 

ball/outer  race  contact 
outer  ring 
spinning  or  sliding 
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In  order  to  understand  the  sudden  failure  phenomenon  in 
angular  contact  ball  bearings,  a time-simulation  analysis 
which  solves  the  equations  of  motion  for  each  ball  is 
performed.  The  analysis  is  performed  on  a 60’  contact  angle, 
90  mm  bore  ball  bearing  which  has  failed  (in  a sudden  manner) 
at  9,000  rpm  during  a test  run.  The  calculation  shows  that  the 
heat  generation  rate  in  this  bearing  increases  significantly 
when  the  bearing  speed  exceeds  a certain  threshold  value.  This 
sudden  increase  of  heat  generation  rate  results  from  ball 
skidding  caused  by  the  gyroscopic  moment  effect.  It  was  also 
found  that  a radial  load  on  the  bearing  can  induce  an  even 
higher  heat  generation  rate,  and  this  radial  load  effect  can 
be  solved  by  the  time-simulation  method.  For  comparison 
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reasons,  the  heat  generation  rate  in  a 15°  contact  angle  ball 
bearing  with  the  same  dimensions  is  also  calculated. 

It  has  also  been  identified  that  the  onset  of  the  ball 
skidding  is  dependent  on  the  friction  in  ball-race  contacts, 
which  are  usually  lubricated  by  elastohydrodynamic  (EHD) 
f . The  traction  coefficients  of  EHD  films  are  dependent 
on  the  sliding  speed.  Three  different  traction  curves  (curves 
describing  traction-sliding  relationship)  are  used  in  the 
analysis  of  ball  motion  and  heat  generation.  It  was  found  that 
traction  coefficients  in  the  thermal  region  are  important  to 
the  onset  of  ball  skidding,  due  to  ball  spinning  which  occurs 
in  an  angular  contact  ball  bearing. 
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CHAPTER  1 
INTRODUCTION 

Scope  of  the  Problem 

Prediction  of  bearing  speed  limit  is  a problem  of 
practical  significance.  In  addition  to  the  centrifugal 
loading,  various  dynamic  effects  may  also  limit  the  bearing 
speed.  A particular  dynamic  effect  which  must  be  considered 
in  the  application  of  angular  contact  ball  bearings  the  ball 
skidding  caused  by  the  gyroscopic  moment  effect. 

The  gyroscopic  moment  is  generated  because  the  ball  self- 
rotational  axis  and  the  ball  orbiting  axis  in  an  angular 
contact  ball  bearing  are  not  parallel.  Ball  skidding  occurs 
when  the  gyroscopic  moment  exceeds  the  friction  torque 
provided  by  the  ball-race  contacts.  The  so-called  minimum 
preload  condition  as  shown  below  is  developed  from  the  above 
fact  [1,2]. 

Mg  < M D (Q.  + Qg)  (1_1) 

where  Mg  is  the  gyroscopic  moment,  D is  the  ball  diameter,  Q. 
and  Qq  are  loads  on  the  inner  and  outer  race  contacts 
respectively,  and  is  the  traction  coefficient. 
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To  calculate  the  heat  generation  rate  resulting  from  the 
ball  skidding  and  to  understand  the  causes  of  sudden  bearing 
failure,  however,  a computer  analysis  is  required. 
Furthermore,  special  care  has  to  be  taken  in  determining  the 
traction  coefficient  in  Eg. (1-1) , not  only  because  the 
elastohydrodynamic  (EHD)  traction  behavior  is  complicated  but 
also  because  there  is  another  sliding  motion,  ball  spinning, 
occurring  in  the  bearing  at  the  same  time. 

Literature  Review 
Methods  of  Bearing  Analysis 

The  simplest  method  for  calculating  the  ball  motion 
(kinematic  velocities)  and  the  corresponding  heat  generation 
rate  is  to  use  the  race-control  theory,  which,  however,  does 
not  consider  a balance  of  forces.  The  race-control  theory  was 
proposed  by  Jones  [3]. 

Harris  [1,2]  was  the  first  to  present  a computer  method 
which  considers  the  force  actions  in  the  calculation  of 
bearing  kinematics.  Harris^  method  is  now  known  as  the  quasi- 
static method.  After  Harris,  two  large-scale  programs,  CYBEAN 
[4]  and  SHABERTH  [5],  were  developed  based  on  the  quasi-static 
method.  However,  the  quasi-static  method,  which  solves  the 
equilibrium  conditions  for  each  ball  as  it  fixed  at  a certain 
position  in  the  bearing,  is  only  valid  when  the  bearing  is 
axisymmetrically  loaded. 


3 


An  improved  method  for  solving  bearing  dynamics  was 
proposed  by  Walters  [6]  and  Gupta  [7,8].  Their  method  was  to 
solve  the  equations  of  motion  by  time  simulation.  By  using 
the  time-simulation  method,  the  time-cycle  variation  of  the 
loads  and  contact  angles  of  a ball  in  a radially  loaded 
bearing  can  be  calculated. 

Sudden  bearing  failure  caused  by  dynamic  effects  has  been 
the  subject  of  many  bearing  analyses.  Besides  the  gyroscopic 
moment  effect,  some  other  effects  also  can  cause  sudden 
failure  in  bearings  including  radial  bearings.  For  example, 
at  high  speeds,  the  centrifugal  force  may  unload  the  ball- 
inner  race  contacts  and  the  drag  resulting  from  the  lubricant 
in  the  bearing  cavity  can  cause  another  type  of  ball  skidding 
(different  from  that  caused  by  the  gyroscopic  moment  effect) 
[9,10,11].  Cage  instability  is  another  example.  At  high 
speeds,  the  smallest  unbalance  in  the  cage  can  lead  to 
significant  ball-cage  interactions.  Also,  when  the  cage  goes 
into  whirl  motion  (a  self-excited  type  of  vibration) 
[12,13,14],  significant  sliding  can  be  induced.  However,  due 
to  the  difficulty,  only  a few  papers  have  reported 
experimental  observation  on  the  ball  motion  [15,16]. 

Elastohvdrodvnamic  Lubrication  Properties 

In  a rolling  bearing,  the  high  contact  pressure  prohibits 
the  development  of  a convectional  hydrodynamic  film.  However, 
a special  form  of  fluid-film  lubrication,  elastohydrodynamic 
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lubrication,  is  able  to  develop  between  the  balls  and  race  of 
a bearing.  The  existence  of  the  EHD  film  is  attributable  to 
the  elastic  deformation  which  creates  a coherent  channel  for 
the  fluid  flow.  Two  particular  issues  in  EHD  lubrication  are 
of  concern  in  this  research:  (1)  the  effects  of  sliding 
(including  ball  spinning)  on  the  EHD  films  and  (2)  the 
calculation  of  the  EHD  traction  coefficient. 

First,  it  was  found  by  Cheng  and  Sternlight  [17]  that 
sliding  (including  spinning)  has  only  a modest  effect  on  the 
film  thickness.  The  significance  of  this  finding  is  that  the 
isothermal  film  thickness  formulas  are  applicable  in  cases 
where  sliding  is  present.  These  isothermal  film  thickness 
formulas  use  the  inlet  viscosity  to  calculate  the  film 
thickness.  It  has  also  been  suggested  that  although  sliding 
does  not  significantly  affect  the  EHD  film  thickness,  the 
temperature  rise  caused  by  sliding  can  lead  to  lubricant 
deterioration  through  chemical  reactions  [ 19 ] . 

The  second  issue,  the  calculation  of  the  traction 
coefficient,  represents  the  most  complicated  aspect  of  the 
EHD  phenomenon.  Unlike  the  film  thickness,  the  traction 
coefficient  is  sensitive  to  the  sliding  speed.  It  is  now  well 
accepted  that  the  lubricant  rheological  behavior  (determining 
traction  force)  can  be  described  by  the  visco-elastic  model 
and  the  limiting  shear  stress  concept  [20-23].  Experimental 
data  on  the  rheological  parameters,  viscosity,  shear  modulus, 
and  limiting  shear  stress  are  available  for  some  common 
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lubricants  [24,25,26].  Since  the  average  pressure  in  the  ball- 
race  contact  reaches  0.8  GPa,  the  lubricant  traction  behavior 
is  governed  by  the  non-Newtonian  effect. 

Objectives  of  Investigation 

The  first  objective  of  the  current  research  is  to 
understand  the  ball  skidding  phenomenon  (due  to  gyroscopic 
moment  effect) , especially  how  the  ball  skidding  leads  to 
sudden  bearing  failure  and  how  the  ball  spinning  affects  the 
onset  of  ball  skidding.  Although  much  work  on  bearing  dynamics 
has  been  performed,  the  ball  skidding  motion  has  never  been 
clearly  described.  Most  of  the  previous  workers  investigated 
bearings  with  small  contact  angles  (less  than  20'),  where  the 
gyroscopic  moment  effect  is  not  as  visible  as  it  is  on  the 
bearings  with  large  contact  angles.  Furthermore,  some  of  these 
investigations  were  based  on  the  quasi-static  solution,  which, 
as  discussed  earlier,  is  accurate  only  when  the  bearing  is 
axisymmetrically  loaded. 

The  second  objective  of  this  research  is  to  study  the 
relationship  between  the  ball  motion  and  the  EHD  traction 
behavior.  This  study  needs  to  calculate  the  working 
temperature  of  the  lubricant  at  the  contacts  and  also  needs 
to  evaluate  how  the  bearing  performance  responds  to  various 
traction  curves.  The  ball  skidding  is  a problem  where  the 
traction  coefficient  can  strongly  affect  the  ball  motion. 
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Introduction  to  Current  Methods 

This  research  evaluates  heat  generation  rates  in  a 15° 
and  a 60°  contact  angle  ball  bearings  (called  Bearing  A and 
Bearing  B)  to  study  the  ball  skidding  phenomenon.  The 
following  methods  are  used: 

1. The  race-control  theory  is  used  in  the  following 
calculations:  (1)  heat  generation  caused  by  ball  spinning;  (2) 
initial  conditions  for  time  simulation;  and  (3)  threshold 
speeds  of  ball  skidding. 

2.  Both  the  quasi-static  method  and  the  time-simulation 
method  are  performed  to  calculate  heat  generation  caused  by 
ball  skidding.  The  time-simulation  method  provides  more 
accurate  solutions  for  non-axisymmetrically  loaded  bearings. 

3.  A hypothetical  traction  curve  (called  Curve  1)  which 
assumes  the  traction  coefficient  is  almost  independent  of 
sliding  speed  (except  in  the  linear  region  of  the  traction 
curve)  is  used  in  most  of  the  calculations.  Two  other 
hypothetical  traction  curves  (Curves  2 and  3)  which  simulate 
true  EHD  traction  behavior  are  used  to  study  which  part  of  the 
traction  curve  really  affects  the  onset  of  ball  skidding. 

4 . Temperature  rise  in  the  EHD  film  caused  by  ball 
spinning  is  also  presented.  The  temperatures  on  the  solid 
boundaries  of  the  EHD  film  are  determined  by  modeling  the 
solids  as  semi-infinite  planes  subjected  to  a moving  heat 
source  [27,28]. 


CHAPTER  2 

MATHEMATICAL  MODEL  FOR  BEARING  ANALYSIS 


Basic  Ecmations  and  Coordinate  Systems 


The  major  purpose  of  bearing  analysis  in  this  work  is  to 
determine  the  heat  generation  in  bearings.  Such  calculation 
is  accomplished  by  solving  the  equations  of  motion  for  each 
bearing  element  (ball  or  race)  to  determine  the  relative 
motions  between  the  elements. 

Ecmations  of  Motions 

According  to  Newton's  law,  the  equations  of  motion  are 
as  follows: 

for  Cartesian  coordinates  (x,y,z), 


m X = S F 


X 


m y = S Fy 
m z = S F^ 


(2-1) 


for  cylindrical  coordinates  (x,r,0), 


m X = S F 


X 


m r = S F 


r 


(2-2) 


about  a fixed  point  which  rotates  with  an  angular  speed  0 


7 
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about  an  axis  parallel  to  the  x-axis,  the  equation  of  motion 
can  be  written  as 


J CJ 

X 

= S 

J {() 

- J 

0 = S 

My 

(2-3) 

J 

+ J 

II 

M 

that 

for 

a ball, 

J = J = J = J,. 

Two  Simplified  Methods 

Quasi-static  method.  The  quasi-static  model  can  be  used 
when  steady-state  conditions  prevail.  This  method  neglects  the 
inertial  forces  (induced  by  accelerations:  x,  y,  z,  if,  $,  w , 
Wy,  and  and  Coriolis  force  (2  m ir  0)  but  considers  the 
centrifugal  force  and  gyroscopic  moment.  The  equations  used 
in  this  method  are  called  force  equilibrium  equations. 

Race-control  theory.  The  race-control  theory  is  based  on 
a set  of  hypothetical  kinematic  constraints  which  provide  a 
simplified  method  for  solving  for  bearing  kinematic 
velocities. 

Coordinate  Systems  for  Bearing  Analysis 

Ball  bearings,  except  for  some  specially  designed 
bearings,  consist  of  an  inner  race,  an  outer  race,  a number 
of  balls,  and  a cage.  Three  coordinate  systems  are  defined  to 
describe  the  position  and  motion  of  each  of  these  elements. 
These  coordinate  systems  are  referred  to  as  inertial 
coordinate  system,  ball-fixed  coordinate  system,  and  contact 
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coordinate  system,  respectively.  These  three  coordinate 
systems  are  shown  in  Fig. 2-1. 


Figure  2-1  Coordinate  systems  for  bearing  analysis. 

Inertial  coordinate  system.  The  inertial  coordinate 
system  is  used  to  describe  the  position  and  motion  (except  for 
the  self-rotation  of  the  ball)  of  each  bearing  element.  The 
race  and  cage  are  described  by  a set  of  rectangular 
coordinates,  denoted  as  (X,Y,Z);  the  balls  are  described  by 
a set  of  cylindrical  coordinates,  denoted  (X,R,0) . The  origin 
of  this  coordinate  system,  point  O,  is  located  at  the  center 
of  the  non-moving  race.  The  X-axis  is  along  the  spindle  axial 
direction. 
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In  the  inertial  coordinate  system,  the  following 
positions  and  velocities  are  defined: 
for  the  balls, 

(Xjj,R|^,0jj)  : coordinates  of  the  ball  center 

: ball  orbiting  velocity  about  the  reference  point 
for  the  inner  race, 

(X^^, : coordinates  of  the  inner  race  center 

0^^:  inner  race  rotating  speed 

Ball-fixed  coordinate  system.  The  ball-fixed  coordinate 
system,  denoted  as  (x,y,z),  is  attached  to  each  ball  center 
to  describe  the  self-rotational  motion  of  each  ball.  The 
x-axis  of  this  coordinate  system  is  parallel  to  the  X-axis  of 
the  inertial  coordinate  system,  and  the  z-axis  points  in  the 
radial  direction.  The  following  velocities  are  defined  in 
these  coordinate  systems: 

Wy,  w^:  ball  spinning  speeds  in  the  x,  y,  and  z 
directions. 

Contact  coordinate  system.  The  contact  coordinate  system 
is  located  at  the  center  of  each  ball-race  contact  ellipse  and 
is  denoted  as  (x,  y,  z)  . In  this  coordinate  system,  the  x-axis 
is  transverse  to  the  rolling  direction,  the  y-axis  is  along 
the  rolling  direction,  and  the  z-axis  is  normal  to  the 
contact.  This  coordinate  system  will  be  illustrated  later. 
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Variables  to  be  Solved 

There  are  nine  variables  to  be  solved,  six  for  the  balls 
and  three  for  the  inner  race  (assuming  the  outer  race  is 
fixed)  . These  are  and  for  the  inner  race  and 

Rjj,  Wy,  and  for  each  ball. 

Sliding  Motions  in  Ball  Bearings 

Most  of  the  heat  generated  in  a rolling  bearing  is  due 
to  sliding  at  ball-race  contacts.  There  are  four  possible 
types  of  sliding  motion  which  may  occur  in  these  contacts. 
The  patterns  of  these  four  types  of  sliding  are  shown  in 
Fig. 2-2 . 

Ball  Spinning 

In  angular  contact  ball  bearings,  the  extensions  of  the 
two  contact  tangents  (Lines  and  in  Fig. 2-3)  do  not 
intersect  at  the  same  point  in  the  bearing  axis.  This  means 
that  angular  contact  ball  bearings  are  prevented  from  a pure 
rolling  condition.  In  effect,  there  is  a relative  spinning 
motion  (ball  spinning)  between  the  ball  and  race  in  the 
contact  area.  The  ball  spinning  friction  will  be  further 
discussed  later  when  methods  of  calculating  bearing  kinematic 
velocities  are  introduced. 
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b)  sliding  due  to 

gyroscopic  moment 


c)  sliding  due  to  the  contact 
surface  curvature 


d)  gross  sliding 


Figure  2-2  Sliding  line  patterns. 


Sliding  due  to  Gyroscopic  Moment  Effect 

In  an  angular  contact  ball  bearing  the  ball  spinning  axis 
is  not  parallel  to  the  ball  orbiting  axis,  as  shown  in  Fig. 2-4 
(where  a thrust  ball  bearing  is  shown) . A moment,  known  as  the 
gyroscopic  moment,  is  required  for  the  precession  motion  of 
balls  (rotation  of  the  ball  spinning  axis) . The  magnitude  of 
the  required  gyroscopic  moment  is 


M 


g 


= j 0b  ^2 
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Figure  2-3  Race  control  theory  and  the  sliding  lines  in 
ball-race  contacts. 

The  required  moment  can  be  obtained  from  the  friction 
forces  at  the  ball-race  contacts.  If  the  friction  forces  do 
not  provide  a sufficient  moment,  the  ball  will  rotate  about 
its  own  y-axis  and  produce  excessive  sliding.  The  well-known 
minimum  preload  condition  can  be  written  as  follows: 

J < 0.5  Djj  (Q.  + Qq)  (2-4) 

where  u^-  is  the  sliding  speed  in  the  x-direction  (transverse 
to  the  rolling  direction) . 
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Figure  2-4  Gyroscopic  moment  in  an  angular  contact 
ball  bearing. 

Sliding  due  to  the  Curvature  of  Contact  Surface 

As  shown  in  Fig.  2 -2c,  the  contact  area  is  a curved 
surface,  thus  distances  from  the  ball  spinning  axis  to 
individual  points  in  the  surface  are  different.  Therefore, 
pure  rolling  occurs  at  most  at  two  points.  The  middle  section 
slides  in  one  direction,  and  the  outer  sections  slide  in  the 
opposite  direction.  For  pure  radial  bearings,  this  sliding 
friction  represents  the  major  friction  source. 

Gross  Sliding  in  Ball-Race  Contacts 

Normally,  gross  sliding  does  not  occur  in  rolling 
bearings  except  due  to  the  following  reasons:  ball  skidding 
due  to  gyroscopic  moment,  ball  skidding  due  to  lubricant  drag, 
acceleration  or  deceleration  of  bearings,  and  collision 
between  cage  and  ball.  The  pattern  of  gross  sliding  is  shown 
in  Fig.2-2d. 
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Force  Ecfuilibrium  Equations  for  Load  Distribution  Analysis 

The  load  distribution  in  a bearing  can  be  calculated  by 
solving  the  force  equilibrium  equations  (for  steady-state 
loads) , and  an  iterative  method  has  to  be  used  to  solve  these 
equations. 

The  force  equilibrium  equations  for  each  ball  (according 
to  Fig. 2-5)  are  as  follows: 
for  x-direction, 

Q^.  sina^  - Qg  sina^  - m cosa^.  + M Qg  cosa^  = 0 (2-5a) 

for  z-direction, 

cosa^  - Qq  costtg  + M Qi  sina^  “ M Qg  sina^ 

+ Fj.  = 0 (2-5b) 

where  F^  is  the  centrifugal  force  which  can  be  calculated  by 

V (2-6) 

In  Eq. (2-5a  and  b) , the  terms  that  include  the 
coefficient  of  friction  n are  the  reaction  moments  against  the 
gyroscopic  moment.  These  terms  are  included  to  deal  with  the 
effects  of  gyroscopic  moment  on  the  load  distribution.  The 
coefficient  of  friction  can  be  calculated  from  the  following 
equation: 

0.5  M djj  ( Q.  + Qq  )=  Mgy  (2-7) 

where  Mg^  is  the  gyroscopic  moment. 


16 


Assuming  that  the  outer  race  is  fixed,  the  equilibrium 
equations  for  the  inner  race  are  as  follows: 


F = 
X 

s 

Qi 

sina^ 

(2-8a) 

'■y  = 

s 

Qi 

cosa^  cos0|^ 

(2-8b) 

s 

Qi 

cosa^  sin0jj 

(2-8c) 

In  the  above  equations,  the  range  of  summation  is  the 
complete  ball  set.  The  calculation  of  the  ball-race 
interaction  forces  (Q^  and  Q^)  requires  a complex  contact 
geometry  analysis  and  is  discussed  in  Chapter  3 . 


Figure  2-5  Forces  and  moments  acting  on  a ball  and  the  races 
(a  simplified  model) . 
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Solving  Bearing  Kinematics  bv  Using  the  Race-Control  Theory 
Race-Control  Theory 

Due  to  the  fact  that  the  ball  self-rotational  axis  is  not 
parallel  to  the  bearing  axis,  the  kinematics  of  angular 
contact  ball  bearings  is  somewhat  more  difficult  than  that  of 
pure  radial  bearings.  The  simplest  method  to  calculate  the 
kinematics  of  this  type  of  bearing  is  to  use  the  race-control 
theory. 

The  race-control  theory,  which  is  developed  from  the 
concept  of  dry  friction,  assumes  that  the  ball  spinning  only 
occurs  on  one  of  the  races.  The  race  on  which  no  spinning 
occurs  is  called  the  controlling  race  (being  the  race  that 
produces  more  resistance  to  spinning) . Usually,  the 
controlling  race  is  the  inner  race  at  low  speeds  and  the  outer 
race  at  high  speeds.  The  race-control  theory  provides  an  extra 
condition  for  solving  the  bearing  kinematics.  However,  for 
bearings  lubricated  by  EHD  films,  the  race— control  theory 
somewhat  departs  from  the  real  situation  because  the  relative 
spinning  can  occur  on  both  races. 

The  case  shown  in  Fig. 2-3  is  an  outer  race-control  case 
where  the  ball  self-rotational  axis  intersects  Line  at  the 
same  point  in  the  bearing  axis.  The  resultant  velocity  of  ball 
orbiting  velocity,  ball  spinning  velocity,  and  the  outer  race 
velocity  (which  is  zero)  is  parallel  to  the  contact  surface, 
and  there  is  no  spinning  on  the  outer  race.  For  the  inner 
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race,  the  net  velocity  is  not  parallel  to  the  contact  surface. 
The  sliding  lines  caused  by  the  relative  spinning  is  also 
shovm  in  the  inner  race  contact  while  the  sliding  lines  in  the 
outer  race  represent  the  sliding  due  to  the  curved  surface. 

By  utilizing  the  race-control  theory,  three  equations  can 
be  written  according  to  Fig. 2-6.  In  the  following  equations, 
Eq. (2-9a)  is  written  according  to  the  no  spinning  condition 
in  the  controlling  race;  Eq.  (2-9b)  is  written  for  the  no  gross 
sliding  condition  in  the  controlling  race;  Eq. (2-9c)  is  for 
the  no  gross  sliding  condition  in  the  non-controlling  race. 
In  these  equations,  subscript  c refers  to  the  controlling 
race,  and  subscript  n refers  to  the  non-controlling  race, 

sina^  + cosa^  + sina^ 

= 0^.^  sina^  (2-9a) 

0.5  d|^  cosa^  - 0.5  dj^  sina^  + (rp  + 0.5  dj^  cosa^)  <p^ 

= (rp  + 0.5  djj  cosor^)  0^^  (2-9b) 

0.5  djj  cosa^^  £J^  - 0.5  dp  sina,^  + (rp  + 0.5  dp  cosa^)  0p 

= (rp  + 0.5  dp  cosa^)  0^^  (2-9c) 

where  rp  is  the  radius  of  the  pitch  circle. 

Calculation  of  Ball  Spinning  Friction 

The  following  equations  can  be  used  to  determine  the  ball 
spinning  velocities  and  rolling  velocities.  The  ball  spinning 
velocities  for  the  inner  race  and  outer  race  are 
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R 


Figure  2-6  Rotational  velocities  in  a ball  bearing. 


"si  " "x  + “z  sina^  - 0.^  sina.  (2-lOa) 

"so  = "x  "z  (2-lOb) 

The  rolling  velocities  are 

"ri  ""  "x  “ “z  <^b  “ <^ir  cosa.  (2-lla) 

"ro  " "x  cosCq  - sina^  + <f>^  cosa^  (2-llb) 

A parameter  which  is  often  used  to  assess  bearing 
performance  is  called  spin-to-roll  ratio.  This  ratio  is 
defined  as  (for  both  inner  race  and  outer  race) 


spin-to-roll  ratio  = 


(2-12) 
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The  frictional  moment  caused  by  the  ball  spinning  motion 
can  be  derived  by  integrating  the  moments  produced  by  friction 
forces  over  the  contact  ellipse.  If  the  pressure  distribution 
in  the  contact  area  is  Hertzian  and  the  coefficient  of 
friction  is  a constant,  this  frictional  moment  will  be 


3^QaE2 


(2-l3a) 


where  E2  is  the  complete  elliptic  integral  of  second  kind  as 
defined  in  Appendix  A,  and  a is  the  length  of  the  semi-major 
axis. 

The  heat  generation  rate  due  to  ball  spinning  can  be 
calculated  by 


H = M 


w . 
SI  SI 


+ M 


(J  • 
SO  SI 


(2-13b) 


Where  subscripts  i and  o refer  to  inner  race  and  outer  race 
respectively. 


Equations  of  Motion  for  a Ball 
The  race-control  theory,  which  does  not  consider  the 
actual  force  equilibrium  conditions,  has  two  drawbacks:  (1) 

the  ball  spinning  may  occur  on  both  races,  and  (2)  the 
gyroscopic  moment  effect  is  not  considered.  This  section 
presents  another  method  which  solves  the  equations  of  motion 
to  determine  bearing  kinematic  velocities. 
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Discretization  of  the  Contact  Area 

Since  sliding  velocity,  pressure,  and  coefficient  of 
friction  vary  from  point  to  point  within  a contact  area,  the 
frictional  force  and  moment  occurring  in  a contact  must  be 
calculated  by  integration  over  the  contact  area.  For  this 
reason,  each  contact  area  is  divided  into  slices  as  shown  in 
Fig.  2-7.  The  following  two  assumptions  are  used  in  the 
discretization; 

l.The  contact  area  is  discretized  in  the  x-direction  only. 
This  is  because  contact  areas  are  often  narrow  ellipses. 

2. Inside  each  slice  the  sliding  velocity,  pressure,  and 
coefficient  of  friction  are  assumed  to  be  constant,  and 
their  values  for  each  slice  are  calculated  at  the  center 
point  of  the  slice. 


y 


u 


sy 


X 


2a 


Figure  2-7  Discretization  of  the  contact  area  for 
the  calculation  of  sliding  friction. 
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Deriving  the  Equations  of  Motion 


To  solve  for  the  ball  motion,  the  following  four  unknowns 
should  be  solved:  and  The  equations  of  motion 
(four  equations  for  each  ball)  which  consider  the  effects  of 
the  inertial  forces  and  inertial  torques  can  be  written  by 
making  use  of  Fig. 2-8  as  follows: 


Figure  2-8  Forces  and  moments  acting  on  a ball  and  the 
races. 

forces  in  rolling  direction, 

••  * * r ^syi 

+ 2 mjj  Rjj  = - J q.  dA^  + 

^si 

r '^syo 

J ^*<“30)  % ^0  + ■'drag  + «cage 

^SO 

where  F^^^g  and  represent  forces  exerted  by  the  lubricant 

and  cage. 
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moments  in  the  x-direction. 


Jt,  = 


r ^syi 

J f»(U5,)  q,  r„,  dA,  + 


SI 


syo 


u 


^czo  ^0 


(2-14b) 


so 


moments  in  the  y-direction, 

Jfa  Wy  - Jjj  = - 0.5  djj  [ J 


^sxi 


u 


qi  dA^.  + 


SI 


r , X sxo 


^0  ^CZO  ^0  ^ (2-14c) 


U 


moments  in  the  z -direction. 


SO 


^syi 


• • r I 

"z  + “y  '»'b  = - 1 "(“si)  ‘Ji  >^0X1  + 

U..Z 


SI 


r ^syo 

J ^^^so^  *^o  ^CXO  ^0  (2-14d) 


u 


so 


In  Eq. (2-14a  to  d) , 

^sxi ' ^sxo  ^ sliding  speed  in  the  x-direction  (for  inner 
race  and  outer  race) 


u$yi ' ’^syo  “ sliding  speed  in  the  y-direction  (for  inner 
race  and  outer  race) 


u • = / u^  — • + — 

SI  sxi  “ syi 
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u 


so 


+ 


^czi ' ^czo'  ^cxi ' ^cxo*  distances  from  a point  in  the 

contact  area  to  the  ball  center  (as  indicated  in 
Fig. 2-9) 

q;  the  load  on  each  slice 
Q:  the  total  load  on  the  whole  contact 
Now  the  sliding  velocities  have  to  determined. 


Figure  2-9  Geometric  model  for  calculating 
sliding  velocities. 

Sliding  Velocity  on  Each  Slice 

The  sliding  velocity  at  a point  (x,y,'z)  in  a contact  can 

be  written  as  follows: 

sliding  in  x-direction, 
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^ ^sxo  " 0*5  Wy  djj  (2-15a) 

(this  sliding  velocity  is  caused  by  the  gyroscopic  moment) 
sliding  in  y-direction, 
for  inner  race, 

^syi  “ ^czi  ~ “z  ^czi  “ ^czi  ^czi  (2-i5b) 

for  outer  race, 

^syo  = ■ ‘^x  "^czo  "z  ^czo  " ^ ^czi  (2-l5c) 

The  distances  from  the  point  at  (x,y,z)  to  the  ball 


center  and 

the 

bearing  axis 

can  be  derived  as  follows 

^cxi 

= 

in 

• 

o 

djj  sina^  + x 

cosa^ 

+ ^ sina^ 

^czi 

= 

0.5 

d|j  cosa^  + X 

sina^ 

- z cosa. 

^czi 

= 

■ 

- 0.5  djj  cosa^- 

- X 

sina^  + z cosa^ 

^cxo 

= 

in 

• 

o 

djj  sina^  + x 

cosa^ 

+ z sina^. 

r* 

czo 

= 

0.5 

djj  cosa^  + X 

sina^ 

- z cosa^ 

R _ _ 

czo 

= 

- 

■ 0.5  d|j  cosa. 

- X 

sina^  + z cosa^ 

The  z coordinate  of  this  point  can  be  calculated  by  the 
following  equation; 

^ “ ^deform  “ deform  “ (2-16) 

where  R^jeform  radius  of  curvature  of  the  contact  surface 

(see  Fig. 2-2c) . 

The  contact  surface  is  a curved  surface.  According  to 
Hertzian  solution,  the  radius  of  the  contact  surface  is  the 
harmonic  mean  of  the  curvatures  of  the  two  contacting  bodies 
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harmonic  mean  of  the  curvatures  of  the  two  contacting  bodies 
as 

111 

= ( + ) / 2 (2-17) 

^deform  0.5  dj^  ^ 

where  f is  the  curvature  factor  of  the  race  ( f is  defined  in 
Chapter  3 ) . 

Since  the  contact  area  is  a narrow  ellipse,  the  curvature 
in  the  rolling  direction  can  be  neglected,  and  only  the 
curvature  in  the  transverse  direction  has  to  be  considered. 

Heat  Generation  Rate 

The  heat  generation  rates  in  each  contact  is 

H = J /i(u^)  q u^  dA  (2-18) 

EHD  Traction  Coefficient 
Typical  Traction  Curve 

The  bearing  analyst  needs  to  know  traction  coefficients 
at  various  sliding  speeds  (the  sliding  speed  varies  within  the 
contact  zone) . Since  the  calculation  of  the  EHD  traction 
coefficient  is  somewhat  complicated,  it  is  necessary  to 
establish  the  traction-sliding  relationship  (traction  curve) 
prior  to  the  bearing  analysis.  A typical  traction  curve  as 
shown  in  Fig. 2-10  exhibits  three  different  regions.  The  first 
region  is  the  linear  region  where  the  linear  relationship  of 
Newtonian  flow  applies.  The  second  region  is  the  nonlinear 
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to  the  non-Newtonian  effect.  The  third  region  is  the  thermal 
region  where  the  traction  coefficient  falls  as  the  sliding 
increases . 

Parameters  in  the  Traction  Curve 

The  traction  curve  in  the  Fig. 2-10  can  be  represented  by 
four  parameters:  u^^,  the  sliding  speed  at  which  maximum 
friction  occurs;  the  coefficient  of  friction  at  zero 
sliding;  the  maximum  coefficient  of  friction;  the 
coefficient  of  friction  at  high  sliding  speeds.  In  the  later 
analysis  of  this  paper,  this  traction  curve  as  well  as  these 
parameters  will  be  referred  to. 


2-10  Typical  traction  coefficient  of 
elastohydrodynamic  lubrication. 


Figure 
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Common  Values  of  Coeff icieirt:  of  Friction 

For  bearing  analysis,  the  selection  of  the  traction 
coefficient  can  make  use  of  Table  2-1  which  shows  typical 
values  of  coefficient  of  friction  for  metal-to-metal  contacts 
under  various  lubrication  conditions. 

Table  2-1  Typical  values  of  coefficient  of  friction 

type  of  lubrication  coefficient  of  friction 

lubricated  by  HD  film  0.005-0.02 

lubricated  by  EHD  film  0.02-0.1 

boundary  lubrication  0.08-0.14 

dry 


0.25-0.4 


CHAPTER  3 

BEARING  GEOMETRIC  MODEL  AND  NUMERICAL  PROCEDURE 

In  bearing  analysis,  the  calculation  of  contact  loads, 
contact  angles,  and  contact  areas  requires  a complex  ball-race 
interaction  model.  This  chapter  presents  such  a model  along 
with  the  complete  numerical  procedure  for  bearing  analysis. 

Ball-Race  Interaction 

Based  on  Harris'  model  [1,2],  this  section  develops 
necessary  equations  for  calculating  the  interaction  caused  by 
preloads,  axial  loads,  radial  loads,  and  tilting  moments. 

Basic  Geometry  of  Ancmlar  Contact  Ball  Bearing 

Some  basic  bearing  geometric  parameters  are  introduced 
here.  Fig. 3-la  shows  a load-free  angular  contact  ball  bearing 
and  the  following  parameters: 
a : contact  angle 

f^  dj^:  length  of  the  inner  race  curvature  radius 
fg  d|^:  length  of  the  outer  race  curvature  radius 
The  two  radii  of  curvature  are  slightly  larger  than  the 
ball  radius  to  ensure  point  contacts  between  the  balls  and 
races.  The  contact  angle  is  the  angle  between  the  contact  line 
(line  connecting  the  contact  point  and  the  ball  center)  and 
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the  R-axis.  Furthermore,  in  Fig. 3-la,  b,  and  c.  Points  O',  B, 
and  I are  called  outer  race  curvature  center,  ball  center,  and 
inner  race  curvature  center,  respectively  (Fig. 3-lb  shows  a 
loaded  ball  bearing) . These  three  centers  are  often  referred 
to  in  the  later  discussion.  Also  shown  in  Fig.  3-la,  d^^  is  the 
distance  between  the  two  race  curvature  centers  under  the 
load-free  condition.  This  distance  is 


^io  ^0 


I:  inner  race  curvature 
center 

B:  ball  center 


c)  relative  positions  of 
I and  B to  outer  race 
curvature  center  0’ 


Figure  3-1  Notations  for  the  calculation  of  contact  angles 
and  deflections  of  a ball:  a)  an  unloaded  ball;  b) 
a loaded  ball;  and  c)  relative  positions  of  Points 
I and  B to  O'. 


Displacement  of  the  Inner  Race 

Fig. 3-2  shows  how  a bearing  inner  race  (the  big  circle) 
responds  to  a radial  load  and  a tilting  moment.  The  following 
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discussion  is  to  determine  the  displacement  of  the  inner  race 
- in  terms  of  the  displacement  of  the  curvature  center  I of 
a ball  located  at  the  position  <p^. 


C^orc  + djocosa,)  cos9’b  sin  0^ 


a) 


b) 


Figure  3-2  Displacement  of  the  inner  race  (the  big  circle) 
due  to  a)  a translation  in  the  Y-Z  plane;  b)  a 
misalignment  angle  (-0^)  about  the  z-axis. 

Some  more  parameters  have  to  defined,  and  they  are 
(X|j,Rjj,0|j)  ; the  coordinates  of  the  ball  center  B. 

(0,  Rqpj./  0[j)  ’ the  coordinates  of  the  outer  race  curvature 
center  O ' , and  where 


^orc  = ■ (fo  " 


the  displacement  of  inner  race  due  to  axial  preload. 


32 


5p^:  the  displacement  of  inner  race  due  to  radial  preload. 
(X^^, : the  new  position  of  the  inner  race  center  when 
the  inner  race  is  siabject  to  a radial  load. 

Therefore,  the  displacement  of  the  inner  race  center  I, 
as  represented  by  51^  and  Sl^,  under  all  possible  loads  can  be 
rewritten  as 

“x  = ■'px  + + <"orc  + <*io  °°=“0> 

■sir  = («pr  + ’'ir  + ^ir  °°=®ni 

If  there  is  no  tilting  moment,  the  tilting  angle  is 

zero  and  the  above  equations  become 
~ ^px  ^ir 

51^  = 5p^  + COS0J,  + sin0jj 

Deformation  Force  in  the  Ball-Race  Contact 

After  the  displacements  of  the  inner  race  center  I,  51^ 
and  51^,  are  determined,  the  equations  for  calculating  the 
interactions  between  a ball  at  position  (X|^,  R|^,  (j>^)  and  the 
two  races  are  developed  here.  Again,  some  new  notation  is  used 
(referring  to  Fig.3-lc); 


'‘ib  “ '^io 

sinttQ  + 

«^x  - 

- ^orc) 

^ib  ” ^io 
“‘bo  = “‘b 

COSUq  + 

<SI,  - 

^orc^ 

“‘bo  = “b 

^orc 

The  deflections  in  the  ball/race  contact  and  the  contact 
angles  then  can  be  calculated  as  follows; 
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For  inner  race. 


tan’^ 

(3- 

-la) 

= 

(f^  - 0.5)  djj 

(3- 

-lb) 

For  outer 

race, 

“o  = 

tan"^ 

(=‘bo/''bo) 

(3- 

-ic) 

^0  = 

cr 

o 

ro 

- - 

(fo  - 0-5)  db 

(3- 

-Id) 

The  elastic  deformation  forces  corresponding  to  the  above 
deflections  are 


Q = K 5^/^  (3-2) 

The  ball-race  contact  stiffness  K is  dependent  on  the 
conformity  of  the  two  contacting  bodies,  and  the  formula  for 
calculating  this  stiffness  is  discussed  in  Appendix  A. 

Numerical  Procedure  for  Load  Distribution  Analysis 

The  load  on  each  ball  in  a bearing  can  be  obtained  by 
solving  the  force  equilibrium  equations  of  Eq. (2-5)  and  (2-8) . 
The  variables  to  be  solved  are  Z^.^,  X|j,  and  R^j.  The 
force  equilibrium  equations  are  non-linear  algebraic 
equations,  therefore,  these  equations  are  solved  by  Newton- 
Raphson's  iteration  method  in  this  paper. 

Also  in  the  calculation,  the  balls  are  assumed  to  have 
no  effect  on  each  other  so  that  the  computation  effort  can  be 
significantly  reduced. 
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Equilibrium  Ecmations  for  Paired  Angular  Contact  Ball  Bearings 
Since  angular  contact  ball  bearings  are  usually  mounted 
in  pairs,  the  force  eguilibrium  equation  of  Eg. (2-8a)  has  to 
include  the  load  caused  by  the  other  bearing.  The  method  used 
here  is  to  add  the  axial  force  caused  by  the  other  bearing  to 
the  external  force  term  in  Eg.  (2-8a)  . This  method  can  be 
expressed  as 

F = F ' + F 
X X xp 

where  is  the  external  force  applied  to  the  bearing, 

Fj^p  is  the  axial  force  caused  by  the  other  bearing. 

For  a pair  of  angular  contact  ball  bearings  preloaded 
against  each  other,  the  centrifugal  force  will  increase  the 
preload  between  them. 

Procedure  of  Quasi-Static  Load  Distribution  Analysis; 

Step  1:  Input  bearing  geometric  data,  material  properties 

(elastic  modulus,  density) , bearing  preload,  radial 
load,  and  bearing  speed. 

Step  2 : Calculate  the  centrifugal  force  by  using  the  free 
contact  angle,  and  also  assume  no  gyroscopic  moment 
effect  (assume  the  coefficient  of  friction  is  zero) . 
Step  3:  Let  and  equal  zero,  . 

Step  4:  For  each  ball,  assume  initial  values  for  Xjj  and  Rjj  and 
compute  and  in  Eg. (2 -5a)  and  (2 -5b) . 
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If  Eq. (2-5a)  and  (2-5b)  are  satisfied  for  each  ball, 
go  to  the  next  step. 

If  either  Eq. (2-5a)  or  (2-5b)  is  not  satisfied,  using 
Newton-Raphson's  method  to  obtain  new  and  R|^  and 
repeat  this  step. 

Step  5:  Take  into  account  the  gyroscopic  moment  effect.  Use 

calculated  and  for  each  ball  to  determine  the 

coefficient  of  friction  by  Eq. (2-8c) . Repeat  Step  4 

once  again  with  this  coefficient  of  friction  value. 

Step  6:  Calculate  resultant  bearing  forces  F , F , and  F by 

X y L 

Eq. (2-8a) , (2-8b)  and  (2-8c)  (neglecting  the  friction 
forces  in  these  equations) . 

If  and  are  zero,  the  value  of  F^  is 

stored  as  F to  account  for  the  force  caused  by  the 
xp 

other  bearing. 

If  Eq. (2-8a  to  c)  are  satisfied,  go  to  next  step. 

If  any  of  Eq.  (2-8a  to  c)  are  not  satisfied,  use 
Newton-Raphson's  method  to  obtain  new  X^^, 
and  repeat  from  Step  4. 

Derivatives  for  the  Newton-Raohson^s  Method 

In  solving  the  equilibrium  conditions  by  the 
Newton-Raphson' s method,  the  derivatives  of  contact  loads  and 
contact  angles  relative  to  the  ball  position  are  needed.  In 
the  following,  the  derivatives  for  the  inner  race  contact  are 
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presented.  The  derivatives  for  the  outer  race  contact  are 
similar. 

Q,/  X„  = - 1.5  5,V2  + r,^2,l/2 

Q,/  R;,  = - 1.5  ^ r,^2)V2 

“i/  =^b  = =‘il/<=‘lb^  + '^ib^' 

V "b  “ - + “"ib^) 

Procedure  for  the  Complete  Bearing  Dynamic  Analysis 
Time  Simulation  Procedure 

The  time-simulation  method  solves  the  equations  of  motion 
by  continuously  integrating  these  equations  over  a small  time 
step.  Time  simulation  often  starts  from  the  initial  conditions 
provided  by  the  race-control  theory  solution.  The  integration 
can  be  performed  by  the  Runge-Kutta  fourth  order  formulas. 

A flow  chart  of  the  complete  bearing  analysis  ( including 
load  distribution  analysis  and  dynamic  analysis)  is  presented 
in  the  next  page. 

Quasi-Static  Analysis 

The  quasi-static  method  neglects  the  inertial  force  terms 
- mjj  R|j  Jjj  J|j  Wy,  and  Jj^  - in  the  equations  of 

motion,  Eq. (2 -14a  to  d)  . The  following  iterative  procedure  is 
used  for  this  analysis: 

1.  Fix  the  loads  and  contact  angles  of  each  ball. 
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2 . Perform  a time  simulation  on  each  ball  until  a steady  state 
is  reached. 

The  justification  of  this  procedure  is  that  if  the  loads 
and  contact  angles  of  a ball  are  constant,  the  time  simulation 
will  eventually  reach  a steady  state. 


Figure  3-3  The  flow  chart  of  the  time-simulation  analysis 
(including  the  load  distribution  analysis  and  kinematic 
analysis  by  the  race-control  theory. 
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no 


Input  bearing  data:  geometry, 
material  properties,  loads,  speed. 


Load  distribution  analysis  for 
obtaining  loads  and  contact  angles 
of  each  ball. 


Use  race-control  theory.  Eg. (2-9a, 
b,  and  c) , to  obtain  ball  kinematic 
velocities  (u^,  and  <p^)  . 


Use  the  results  of  the  race-control 
theory  as  initial  conditions  to 
start  the  simulation. 


Integrate  the  equations  of  motion 
Eg. (2-14a  to  2-14d)  over  a small 
time  step  to  obtain  new  ball 
position  and  kinematic  velocities. 


yes 


Compute  heat  generation  rates  from 
the  kinematic  velocities. 


( ) 


CHAPTER  4 

TEMPERATURE  DISTRIBUTION  AND  TRACTION  WITHIN 
ELASTOHYDRODYNAMIC  FILMS 

As  discussed  earlier,  the  traction  coefficient  between 
balls  and  races  has  a significant  influence  on  the  ball 
motion.  This  chapter  discusses  the  procedure  for  calculating 
the  traction  coefficient  of  EHD  films.  There  are  three 
sections  in  this  chapter.  Since  lubricant  viscosity  is 
dependent  on  temperature,  equations  for  a fast  computation  of 
the  temperature  in  the  EHD  film  are  derived  in  the  first 
section.  The  film  thickness  formulas  and  the  effects  of 
sliding  on  the  EHD  film  are  discussed  in  the  second  section. 
The  third  section  presents  lubricant  rheological  models  and 
a flow  chart  for  the  calculation  of  traction  coefficient. 

Bearing  Bulk  Temperature  and  Temperature  in  the  EHD  Film 
Bearing  Bulk  Temperature 

The  most  important  aspect  of  the  bearing  bulk  temperature 
is  that  it  determines  the  EHD  inlet  temperature  and 
consequently  the  film  thickness.  However,  the  calculation  of 
the  bearing  bulk  temperature  is  difficult  because  the  heat 
transfer  coefficients  required  for  such  calculation  are 
usually  unknown.  The  following  equation  is  written  for  an 
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estimation  of  this  temperature.  This  equation  assumes  that  all 
the  heat  generated  in  a bearing  is  taken  away  by  a lubricant 
flow. 

H 

"^bulk  “ "^oil  (4-1) 

^oil  ^oil  ® 

where  H is  the  heat  generation  rate  of  the  bearing,  and  Q is 
the  oil  flow  rate. 

Energy  Ecmation 

Temperature  distribution  can  be  calculated  from  the 
energy  equation  as  follows: 

d^e  du  2 dd 

k + r]  ( ) = p C (4-2) 

dz^  dz  3t 

In  this  equation,  the  two  terms  on  the  left-hand-side 
represent  the  heat  conduction  in  the  z-direction  and  the  heat 
generated  due  to  shearing  of  the  fluid  in  the  film 
respectively.  The  right-hand-side  term  represents  the  change 
in  the  fluid  internal  energy. 

The  specific  meaning  of  the  energy  equation  is  that  the 
heat  generated  in  the  contact  is  either  conducted  to  the  metal 
boundaries  or  convected  away  by  the  fluid. 

Solid  Surface  Temperature 

The  temperature  on  the  solid  boundary  surface  can  be 
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calculated  by  the  following  equation  developed  by  Archard 
[27]: 

de 

- k ( ) 

1 y dz  z=0  or  h 

^(y)  = - y f , (4-3) 

y 7T  V -b  J Y - 


where  subscript  s represents  the  solid. 

This  equation  can  be  derived  by  the  analogy  of  a heat 
source  moving  along  a semi-infinite  surface  [13].  Note  that 
V is  the  rolling  speed. 

It  also  can  be  assumed  that  the  heat  transferred  into  the 
solid  surface  is  at  a constant  rate  over  the  contact  area; 
then  the  integration  in  Eg. (4-3)  can  be  carried  out  to  give 

ae 

k ( ) y y + b 

dz  z=0  or  h 

«s(y)  - jr-'  — (4-4) 

y IT  />5  Cj  kj  V 


Temperature  Profile  Obtained  from  the  Energy  Equation 

If  shearing  takes  place  in  the  fluid  film,  the 
temperature  increase  in  the  film  can  be  determined  from  the 
energy  equation  as 


(1+E)  d = + c,  z + c- 

2 k ^ ^ 

where  q is  the  specific  heat  generation  rate,  and 

2 

p c V z 

E = (4-5) 

2 k Ay 
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E is  a factor  to  calculate  the  heat  convected  away. 
For  the  boundary  conditions 

^ at  z = 0 and  z = h 

we  find 


q 

B = z (h-z)  + 

2 k (1+E) 


(4-6) 


dp/dy  = 0 


ball 


Figure  4-1  Distributions  of  the  velocity,  shear  rate, 
stress,  and  temperature  across  an  EHD  film  due 
to  sliding. 


Shear  Stress  Distribution  rslidina  Occurs) 

To  calculate  the  heat  generation  rate,  the  shear  stress 
distribution  across  the  film  should  be  determined. 

dr  dp 

dz  dy 

and  according  to  the  boundary  condition  (Fig. 4-1) 


r = r] 


h 


, at  z = 


2 
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Since  the  shear  stress  varies  linearly  along  the  z- 
direction,  the  shear  stress  due  to  the  pressure  gradient  is 


dp  h 

T = ( z)  + rj 

dy  2 


^s2  - ^sl 


Heat  Generation  due  to  Sliding 

In  the  Hertzian  contact  zone,  the  heat  generation  due  to 
the  pressure  gradient  is  small.  It  is  reasonable  to  assume 
that  all  the  heat  is  generated  from  the  sliding.  In  such 
cases,  the  strength  of  the  heat  source  across  the  film  is 
constant. 

If  0^  represents  the  temperature  on  the  surface  and 
is  the  temperature  at  the  center  of  the  film,  the  following 
relationship  can  be  written, 
dv 

q = r = constant 

dz 

Thus,  when  z = h/2, 

1 3v 

0 = T ( ) z (h-z)  + 0 

2k  (1+E)  dz 

30  1 3v 

= r ( ) (h-2z) 

dz  2k  (1+E)  dz 


for  h = h/2,  and  the  temperature  at  the  central  plane,  0^  can  be 
obtained. 

1 dv  h^ 

0c  = T ( ) + 0^  (4-7) 

2k  (1+E)  dz  4 
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and  now  E=  (pCVh^/8k  Ay) 

For  the  calculation  of  the  heat  transferred  to  the  solid 
boundary  is 

96  1 dv 

- k I 1 = T ( ) h 

dz  z=0  or  h 2 (1+E)  dz 

Temperature  Increase  due  to  Pressure  Gradient  (for  no  Sliding) 
Sliding  represents  the  major  heat  source  in  a bearing  and 
the  temperature  rise  caused  by  sliding  is  much  higher  than 
that  caused  by  pure  rolling.  When  there  is  no  sliding,  heat 
generation  is  due  to  the  pressure  gradient  and  is  only 
significant  in  the  inlet  region.  The  temperature  distribution 
for  a no-sliding  case  is  presented  here. 

By  following  the  similar  procedure  in  the  previous 
section  and  assuming 

9=6^  at  z = 0 and  z = h 
then  the  temperature  distribution  becomes 

q 

6 = (h-z)  z + ^ 

2 k (1+E) 

where  q is  the  specific  heat  generation  which  can  be 
determined  by 

dp  h dv 

T = ( Z)  

dy  2 dz 


since 
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q = T 


we  can  write 


3v 


dz 


[- 


dp 


( z)  ] 


2 k r?  dy 


By  substituting  the  above  heat  generation  rate  into  the 
temperature  distribution,  the  following  relationship  can  be 
written, 


9 = 


dp  h 2 

■[  ( z)  ] z (h-z)  + 9^ 


2 k r?  (1+E)  dy  2 


In  order  to  obtain  the  location  where  the  maximum 
temperature  is,  the  temperature  is  differentiated  over  z. 
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dz 


dp  2 

( ) [(h-2z)(8z^  - 8zh  + h^) ] 


(2k)4r?  (1+E)  dy 


Thus,  the  maximum  temperature  occurs  at  the  following 
positions  (as  shown  in  Fig. 4-2); 

2 ± J~2 


z = 


dp  2 1 


max 


( ) 


(4-8) 


(2k)  477  (1+E)  dy  16 


and  now  E = ( (3-27?)  />  C V/  16  k Ay) 
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For  the  calculation  of  the  heat  transferred  to  the  solid 
boundary  is 


-k 


dff 


dz  z=0  or  h 


1 dp  2 

( ) h 

(2k)4r7  (1+E)  dy 


for  zero  sliding  but 
dp/dy  ^ 0 

ball 


race 


Figure  4-2  Distributions  of  the  velocity,  stress,  and 
temperature  across  an  EHD  film  due  to  a 
pressure  gradient  (sliding  is  zero) . 

EHD  Film  Thickness  and  Effects  of  Sliding 
This  section  presents  the  film  thickness  formulas  and 
discusses  briefly  how  these  formulas  are  derived.  It  has  been 
found  that  sliding  has  only  a modest  effect  on  the  film 
thickness.  An  explanation  of  this  EHD  property  is  also 
presented. 

The  Reynolds  Ecmations 

The  primary  purpose  of  lubrication  is  to  reduce  friction 
and  wear  between  two  rubbing  bodies  by  inserting  a low 
shearing  material  between  them.  The  lubrication  effects  of  an 
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elastohydrodynamic  film  are  accomplished  by  the  fluid  actions 
governed  by  the  Reynolds  equation  (one-dimensional) : 


3 3p  3h 

( ) = 12  (Vj  + V2)  (4-9) 

3y  Tj  dy  dy 


The  derivation  of  the  Reynolds  equation  and  the  method 
for  solving  this  equation  are  presented  in  Appendix  B. 


Film  Thickness  for  EHD  Lubrication 

The  most  important  parameter,  the  film  thickness,  can  be 
calculated  by  using  the  formulas  presented  by  Hamrock  and 
Dowson.  The  following  film  thickness  formulas  are  known  as  the 
isothermal  film  thickness  formulas  because  they  are  obtained 
by  assuming  the  viscosity  is  constant  across  the  z-direction. 
The  minimum  film  thickness 


= 3.63 


uO.68  qO.49  g-0.68/Cj 


R 


W 


0.063 


(4-10) 


The  central  film  thickness, 

h.. 


= 2.69 


uO-67  qO.53 


W 


0.067 


(4-11) 


where  G,  R,  U,  W and  k are  non-dimensional  parameters 
G;  elastic  parameter 

R:  equivalent  radius  (defined  in  Appendix  A) 

U : speed  parameter 

W:  load  parameter 
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k:  ratio  of  contact  area  axes 

and, 


U 


r?o  (uj  +U2) 

2 E'  R 


W = 


Q 


E' 


R 


2 


G = a E' 

where  t]q  : the  inlet  viscosity 

a : the  viscosity-pressure  factor  (generally, 

. o 

IS  between  0.01  and  0.02  mm  /N) 


its  value 


Derivation  of  the  Film  Thickness  Formulas 

The  film  thickness  formulas  are  derived  by  solving  the 
Reynolds  equation,  elasticity  theory,  and  lubricant  rheology 
for  various  sets  of  values  of  the  EHD  parameters  (G,R,U,W,/c)  . 
The  exponent  of  each  parameter  is  determined  by  the  least- 
square  curve  fitting  for  exponential  function. 

Film  Thickness  for  Hvdrodvnamic  fHD)  Lubrication 

It  is  very  difficult  to  solve  the  Reynolds  equation,  but 
for  a simpler  case,  the  HD  lubrication,  the  film  thickness  can 
be  more  or  less  observed  from  the  integrated  Reynolds 
equation.  The  integrated  Reynolds  equation: 


3p 


h - 


ay 


= 12  rj  V 


(4-12) 
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where  is  the  film  thickness  at  a location  where  maximum 
pressure  occurs  (dp/dy  = 0) . 

It  can  be  found  that  the  above  equation  is  similar  to  the 
film  thickness  formula  for  HD  Ivibrication: 

h 2 7?  V 

( ) = C (4-13) 

R P 

where  C is  a constant. 

By  comparing  Eq.  (4-10)  and  Eq.  (4-13)  , it  can  be  found 
that  the  major  difference  between  the  EHD  and  HD  lubrications 
is  that  load  has  a very  small  effect  on  the  EHD  film 
thickness. 

EHD  Pressure  Profile  and  Film  Shape 

The  key  reason  that  an  EHD  film  can  develop  in  a highly 
loaded  contact  is  that  the  pressure  of  the  fluid  flow  can 
create  a coherent  channel  for  the  fluid  flow  itself.  The  EHD 
lubrication  thus  has  the  following  two  features  (referring  to 
Fig. 4-3) ; 

First,  the  film  should  have  a constant  thickness  over  the 
high  pressure  area,  and  the  pressure  distribution  should  be 
similar  to  the  Hertzian  pressure  distribution.  The  film 
thickness  has  to  be  constant,  i.e.,  the  channel  has  to  be 
flat,  because  the  viscosity  in  the  contact  region  is  high  and 
a small  gradient  of  the  film  shape  can  result  in  a large 
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pressure  gradient.  In  order  to  create  a flat  film  shape,  the 
pressure  distribution  has  to  be  Hertzian. 


Figure  4-3  Pressure  profile  and  film  shape  of  an  EHD  film. 

Second,  there  is  a pressure  spike  near  the  exit  and  a 
constriction  of  the  channel  immediately  following  the  pressure 
spike.  The  pressure  spike  occurs  because  the  flow  in  the  exit 
has  to  satisfy  the  Reynolds  boundary  condition  which 
stipulates  that  sub-ambient  pressure  (the  cause  of  cavitation) 
should  not  occur: 

dp/dy  =0,  at  the  exit 

According  to  this  boundary  condition,  the  pressure 
profile  must  fall  inside  the  Hertzian  profile  in  the  exit 
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region.  A pressure  spike  thus  occurs  to  produce  two  effects: 
first,  it  causes  the  fluid  pressure  to  fall  rapidly  to  meet 
with  the  boundary  condition;  second,  it  makes  up  the  pressure 
difference  between  the  fluid  pressure  profile  and  the  Hertzian 
profile  in  the  exit  region.  Again,  in  order  to  create  a 
pressure  spike,  the  channel  has  to  be  narrowed.  The  film 
thickness  at  the  constriction  is  called  the  minimum  film 
thickness. 

Effects  of  Sliding  on  Film  Thickness  and  Pressure  Soike 

An  important  property  of  the  EHD  film  is  that  sliding 
has  only  a moderate  effect  on  the  film  thickness.  This 
property  can  be  seen  from  the  film  thickness  formulas. 

First,  it  can  be  seen  that  the  viscosity  used  in  the  film 
thickness  calculation  is  the  inlet  viscosity  - tJq.  This  means 
that  although  the  temperature  in  the  contact  region  may  be 
high,  the  film  thickness  will  not  be  affected  by  this 
temperature.  This  EHD  property  can  be  explained  by  the  fact 
that  the  deformation  of  the  metal  surfaces  and  the  fluid 
action  are  mutually  dependent,  thus  change  of  viscosity  in  the 
contact  zone  only  affects  the  slopes  of  the  film  shape  but  not 
the  thickness.  It  can  further  be  inferred  that  an  EHD  film 
rarely  fails  because  the  inlet  temperature  is  commonly  low. 

Second,  the  speed  used  in  the  film  thickness  formulas  is 
called  the  entraining  speed  which  is  the  average  of  the  two 
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surface  speeds.  Since  in  rolling  bearings  the  predominant 
motion  is  rolling,  the  sliding  contribution  is  small. 

As  for  the  pressure  spike,  it  can  also  be  expected  that 
sliding  does  not  have  a significant  influence  on  it.  This  is 
because  the  pressure  spike  is  generated  to  satisfy  the 
Reynolds  boundary  condition  (to  avoid  cavitation) . Sliding 
does  not  change  the  boundary  condition. 

Lubricant  Deterioration 

Because  lubricants  have  relatively  low  conductivity,  high 
temperatures  in  the  lubricating  film  are  expected.  Besides 
affecting  the  working  viscosity,  high  temperature  may  lead  to 
significant  chemical  reactions  (oxidation  or  polymerization) 
and  lubricant  deterioration.  It  has  been  pointed  out  that 
lubricant  deterioration  becomes  critical  when  the  temperature 
exceeds  66 °C  (the  rate  of  oxidation  doubles  for  each  10 'C  rise 
in  temperature  of  the  oil  above  66 °C).  The  chemical  reactions 
can  develop  increased  viscosity,  acids,  carbon  residue,  and 
sludge.  However,  since  lubricant  passes  through  the  contacts 
in  a short  period  of  time,  lubricant  deterioration  is  not  a 
problem  of  concern  unless  only  a trace  of  lubricant  is  used. 

Non-Newtonian  Models  and  Procedure  for  Calculating  EHD 

Traction  Coefficient 

The  EHD  traction  behaves  as  described  by  the  traction 
curve  in  Fig. 2-10.  The  calculation  of  such  traction  curve 
requires  detailed  lubricant  rheological  data.  Since  it  is 
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cumbersome  to  calculate  the  traction  coefficient  while  solving 
the  equations  of  motion  for  a bearing,  the  traction  curve,  as 
suggested  by  Gupta  [3],  is  better  to  be  determined  before 
performing  bearing  analysis. 

Non-Newtonian  Behavior 

The  Newtonian  fluid  model  as  shown  below  is  not 
applicable  in  the  traction  calculation  for  EHD  lubrication: 


The  calculation  of  traction  should  use  non-Newtonian 
models  as  introduced  below.  The  non-Newtonian  behavior  arises 
from  the  fact  that  the  increase  of  lubricant  viscosity  causes 
the  lubricant  to  behave  as  a solid.  A solid  may  display 
elastic  and  plastic  behaviors  which  are  the  causes  of  the  non- 
Newtonian  behavior.  The  following  is  a frequently  used  non- 
Newtonian  model. 

Visco-Elastic  Model  and  Nonlinear  Viscous  Model 

The  combined  solid  and  fluid  (viscous)  behavior  can  be 
described  by  a model  which  is  a combination  of  the  Maxwell 
visco-elastic  model  and  the  nonlinear  viscous  model.  This 
model  states  that  the  total  strain  rate  of  the  fluid  element 
is  the  sum  of  the  elastic  and  viscous  strain  rates. 
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T 


(4-14) 
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tanh 
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(4-15) 


G*  dt 
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where  G*  = fluid  shear  modulus  (N/m  ) 

G*  = lubricant  shear  modulus  = E/(2  (1+j/)) 

Yg  = elastic  stain  rate  (1/sec) 

= viscous  strain  rate  (1/sec) 

T-|  : the  limiting  shear  stress 

The  elastic  deformation  rate  in  Eq.  (4-15)  is  usually  much 
smaller  than  the  viscous  strain  rate  and  can  be  neglected. 
Thus,  Eq. (4-15)  becomes 

"1  -1  " 

7 = tanh  (4-16) 

r\ 

It  can  be  noted  that  the  viscous  term  is  treated  by  a 
non-linear  function  since  the  viscosity  of  lubricant  behaves 
non-linearly . Ree-Eyring  (1955)  has  proposed  a •'  sinh  " 
function  to  describe  the  nonlinear  viscous  behavior.  The  " 
tanh’^  " is  a more  recent  model  proposed  by  Bair  and  Winer 
[20] . 

The  limiting  shear  stress  is  usually  determined  from 
experiments  by  using  the  following  expression. 

= C P - Tq  (4-17) 

where  P is  pressure,  and  C is  a constant  which  has  to  be 
determined  from  experiments.  Eq.  (4-17)  is  a form  that  is 
suitable  for  experiments. 
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Calculation  of  Viscosity 

Viscosity  is  a fluid  property  which  is  defined  as  the 
resistance  of  fluid  to  shear  motion.  Viscosity  is  very 
sensitive  to  pressure  and  temperature.  There  are  two  models 
which  can  be  used  to  determine  the  lubricant  viscosity.  The 
first  model  is  proposed  by  Barus  as  follows: 

^(Pf^)  = »?o  exp[ap+/3(5Q-5)  ] (4-18) 

where  tJq  is  the  dynamic  viscosity  at  atmospheric  pressure,  a 
is  the  viscosity-pressure  factor,  p is  the  viscosity- 
temperature  factor,  and  $q  is  a reference  temperature. 

Due  to  the  fact  that  the  Barus  formula  gives  extremely 
high  viscosity  at  high  pressures,  another  equation  has  been 
proposed  by  Roelands,  who  used  his  physical-chemical 
background  to  obtain  (1966) 

logiQ(»7  ) + 1.2  9q  + 135  Sq  p Z, 

= (-ii ) (1  + ) (4-19) 

logio(r?o)  + 1.2  9q  + 135  2000 

The  non-dimensional  parameters  Zj  and  Sq  can  be  calcu- 
lated from  the  coefficients  a and  p in  the  Barus  formula. 


a 


Zj  = 


5.1  X 10'®  ( ln(r?Q)  + 9.67  ) 


i8  (0Q  +135) 


(4-20) 


( ln(r?Q)  + 9.67  ) 


(4-21) 
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After  some  manipulations,  the  Roelands  viscosity  can  be 
expressed  as 


0 +135 

T]  = t]q  exp{  [ln(rjQ)+9.67]  [-1+  (1+5.1  x 10*^  p)^^  ( )‘^°]} 

5q+135 


(4-22) 


Experimental  Data  for  SAE  30  Motor  Oil 

Viscosities  of  SAE  30  Motor  oil  under  various  conditions 
(measured  by  the  falling  ball  viscometer  [25])  are  listed  in 
Table  4-1.  The  reference  condition  is  chosen  as  40 "C  and  1 
atm,  the  reference  viscosity  is 
r)g  = 0.071  Pa*  sec 

The  values  of  a and  )3  can  be  obtained  by  interpolating  the 
experimental  data  in  relation  to  the  Barus  formula,  and  the 
results  are 

a = 1.5  X 10’®  (1/Pa) 

)9  = 0.008  (l/'C) 

Comparison  Between  the  Barus  Formulas  and  the  Roelands  Formula 
The  parameters  in  the  Roelands  formula  can  also  be 
determined  as 
Zj  = 0.42 

Sq  “ 0.2 

Viscosity  values  under  different  pressures  calculated  by 
Barus  and  Roelands  foraulas  for  the  temperature  at  40"C  are 
listed  as  in  Table  4-2.  The  Roelands  formula  is  an  improved 
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formula  which  is  considered  to  be  more  accurate.  It  can  be 
seen  that  the  results  of  Barus  formula  are  unreasonably  high. 


Table  4-1  Measured  viscosity  data  for  SAE  30  motor  oil  [25] . 


viscosity 
(lO'^  Pa»sec) 

temperature 

rc) 

pressure  (MPa) 

40 

100 

150 

Atm. 

71.0 

8.98 

3.07 

34.5 

144.0 

15.4 

5.35 

69.0 

291.0 

27.2 

8.29 

138.0 

1222.0 

65.4 

17.1 

207.0 

3766.0 

- 

- 

276.0 

10817.0 

325.0 

57.2 

345.0 

- 

- 

- 

414.0 

- 

1144.0 

185.0 

524.0 

- 

- 

366.0 

552.0 

— 

6130. 

— 

Table  4-2  A comparison  between  the  calculated  viscosity  values 
from  the  Barus  formula  and  Roelands  formula. 


pressure  fN/m  1 

0.5  X 10^ 
g 

1 X 10^ 

2 X 10^ 


n (Barus) 


1808 


3.36  X 10' 


1.07  X 10 


13 


^0 

^0 


r?  (Roelands) 
139  T]q 
2947  T]q 
92900  T]q 


Procedure  for  Calculating  Traction  Coefficient 

The  procedure  of  calculating  traction  coefficient  is 


outlined  in  Fig. 4-3. 
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Figure  4-4  The 
coefficients. 


flow  chart  for  calculating  traction 


CHAPTER  5 

BEARINGS  UNDER  INVESTIGATION 

Bearing  Dimensions  and  Catastrophic  Failure 
Bearing  Dimensions 

The  bearings  investigated  in  this  work  are  two  bearings 
used  on  two  different  milling  spindles.  The  major  difference 
between  these  two  bearings  is  in  their  contact  angles  - one 
is  15°,  and  the  other  is  60°.  Throughout  this  paper,  the 
following  designations  are  used; 

Bearing  A:  a 15°  contact  angle  ball  bearing,  90  mm  bore. 
Bearing  B:  a 60°  contact  angle  ball  bearing,  90  mm  bore. 

Bearing  internal  geometry  varies  from  company  to  company, 
thus  for  reasons  of  comparison,  it  is  justified  to  assume  that 
the  two  bearings  have  the  same  internal  geometry  except  for 
the  contact  angle  (in  fact,  contact  angle  is  the  most 
important  factor  affecting  heat  generation)  . The  following  are 
the  data  for  these  bearings  which  data  are  typical  for 
bearings  of  7018  series  (90mm  bore) . 

ban  diameter  - 12  mm 

inner  race  curvature  factor  f^  = 0.52 
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outer  race  curvature  factor  = 0.52 
pitch  diameter  dp  = 115  mm 
bore  diameter  = 90  mm 
no.  of  ball  = 20 

ball  density  = 7800  kg/m^  (steel) 

ball  mass  m|j  = 0.007  kg 

. . -7  ? 

ball  inertia  = 1.02  x 10  kg-m 

c p 

bearing  material (steel)  E = 2.08  x 10  N/mm 
Poisson's  ratio  = 0.3 

Catastrophic  Failure  of  Bearing  B 

When  the  speed  of  Bearing  B reached  9,300  rpm  in  a test 
run,  it  failed  suddenly;  below  that  speed,  this  bearing  ran 
stably.  This  bearing  is  a double  row  bearing,  and  the  failure 
only  occurred  in  one  row  - the  less  lubricated  row  (this 
bearing  is  mounted  on  a vertical  spindle,  thus  more  flow  goes 
into  the  lower  row  than  into  the  upper  row).  Fig. 5-la  is  a 
photograph  of  the  undamaged  row,  and  Fig. 5-lb  shows  the 
damaged  row.  It  can  be  seen  that  craters  were  formed  on  some 
of  the  balls  in  the  failed  row.  This  indicates  that  the 
bearing  temperature  reached  a very  high  level. 

Bearing-Spindle  Systems  and  Bearing  Bulk  Temperature 
This  section  shows  the  configurations  of  the  spindles 
which  contain  the  two  bearings  to  be  investigated.  These  two 
spindles  have  been  run  continuously  at  6,000  rpm  for  at 
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b) 


Figure  5-1  Photographs  of  the  angular  contact  ball  bearing  in 
Spindle  1:  a)  the  undamaged  row  (due  to  better  lubri- 
cation) ; b)  the  damaged  row,  craters  on  ball  can  be  seen 
on  the  top  several  balls.  Failure  occurred  when  this 
bearing  was  accelerated  to  9,300  rpm. 
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least  an  hour.  The  operating  temperatures  on  the  outer  races 
of  these  two  bearings  have  been  recorded. 

Bearing  A and  Spindle  1 

This  spindle  consists  of  five  15°  contact  angle  ball 
bearings  (Bearing  A represents  the  three  bearings  on  the  tool 
end)  as  shown  in  Fig. 5-2.  Small  contact  bearings  allow  the 
spindle  to  operate  at  high  speeds.  Since  15°  contact  angle 
bearings  have  relatively  low  axial  stiffness,  a large  number 
of  bearings  (5  bearings)  are  used. 

After  running  at  6,000  rpm  for  one  hour,  the  temperature 
rise  on  the  outer  race  was  10 °C.  The  bearings  were  lubricated 
by  grease. 

Bearing  B and  Spindle  2 

Bearing  B (60°  contact  angle)  is  used  in  Spindle  2 to 
support  axial  loads.  This  spindle  is  equipped  with  two 
cylindrical  roller  bearings  which  provide  substantial  radial 
stiffness  (see  Fig. 5-3). 

After  running  at  6,000  rpm  for  an  hour,  the  temperature 
rise  on  the  outer  race  was  20 °C.  This  bearing  was  lubricated 
by  circulating  oil  flow. 

Calculation  of  the  Hertzian  Parameters 

The  equivalent  radii  and  Hertzian  parameters  (a*,  b*, 
and  5*)  which  are  needed  for  the  load  distribution  analysis 
are  calculated  in  this  section  using  the  prescribed  bearing 
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Figure  5-2  Spindle  1,  consisting  of  five  15  ‘ contact  angle 
ball  bearings. 


Figure  5-3  Spindle  2,  consisting  of  two  cylindrical  roller 
bearings  and  a 60*  contact  angle,  double  row  ball 
bearing. 
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geometry  data  and  related  equations  in  Appendix  A. 
section  only  shows  the  calculation  for  Bearing  A. 


Calculation  of  the  Hertzian  Parameters 
For  inner  race, 


rb^  = 0.5  Db 

= 6 mm 

f 
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If) 
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II 

= 6 mm 
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= - 6. 
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Making  use  of  the  Eq. (A-17a  to  c) , we  obtain, 
K = 1.0339  ( — — )0-6360  ^ 8.8051 

Ej  = 1.5277  + 0.6023  ln(R  / R ) = 3.556 


0.5968 


E2  = 1.0003  + 


= 1.0209 


(Rx/^y) 


2 

2 E 
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b*  = ( 


■ ) = 3.6935 
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IT 


= 0.6129 


( 


7T 


for  outer  race 


rjjj^  = 0.5  Djj  = 6 mm 


r|jy  = 0 . 5 Djj  = 6 mm 


“6.24  mm ^ ^ry  — 63.3  mm 

thus , 

Rj(  = 156  mm  , Ry  = 6.625  mm  , R = 5.170 
Similarly,  the  following  Hertzian  parameters  can  be  obtained: 
K = 7.7096  , Ej  = 3.4304  , E2  = 1.0256 
a*  = 3.3856  , b*  = 0.4391  , S*  = 0.6450 

Including  the  Material  Data  into  Hertzian  Formulas 

The  material  properties  can  be  incorporated  into  the 
Hertzian  formulas,  Eq. (A-8) , (A-9) , and  (A-10) . If  the 

material  is  steel, 

E=  2.08  X 10^^  N/m^,  U = 0.3 


and 


— = 2.28  X 10^^  N/m^ 

1 - 


The  Hertzian  formulas  become 


= 2.14  X 10^  r1/2 

a = 0.0236  a*  (Q  R) 
b = 0.0236  b*  (Q  R) 


CHAPTER  6 

RESULTS  AND  DISCUSSION 

It  should  be  recalled  here  that  the  two  bearings  under 
investigation  are  (discussed  in  Chapter  5) 

Bearing  A;  free  contact  angle  Cq  = 15° 

Bearing  B:  free  contact  angle  Cq  = 60° 

The  friction  in  the  ball-race  contact  is  assumed  to  obey 
the  traction  curve  as  illustrated  in  Fig. 2-3.  Based  on 
available  experimental  data  [10,11,12],  the  parameters  in  the 
traction  curve  are  assumed  to  be  (two  other  traction  curves 
will  be  assumed  later) 

Traction  Curve  1: 

Ujjj  = 0.2  m/sec 

Mo  = 0 

M„,  = 0.054 
= 0.054 

It  can  be  seen  that  this  traction  curve  assumes  traction 
coefficient  levels  off  at  high  sliding  speeds. 

The  analysis  includes  four  steps.  First,  the  loads  and 
contact  angles  of  each  ball  are  calculated.  Second,  the  race- 
control  theory  is  used  to  compute  the  ball  spinning  friction 
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and  the  heat  generation  over  a range  of  the  speed.  Third,  the 
time-simulation  method  is  used  to  investigate  the  ball 
skidding  phenomenon  caused  by  the  gyroscopic  moment  effect. 
Lastly,  the  lubrication  film  thickness  and  temperature 
distribution  are  calculated. 

Load  Distribution 

The  loads,  the  contact  angles,  and  the  contact  areas  of 
each  ball  have  to  be  determined  by  solving  the  force 
equilibrium  equations,  Eq. (2-4) . This  section  shows  calculated 
results  for  one  particular  case  for  each  bearing. 

For  Bearing  A at  10.000  RPM 

For  the  prescribed  loading  conditions,  the  calculated 
results  are  listed  as  follows: 
axial  preload  : 1,000  N 

bearing  speed  : 10,000  rpm 
operating  contact  angles; 

= 21.06”  , ttg  = 13.96” 

loads: 

= 163.7  N , Qg  = 251.6  N 
contact  areas; 
for  inner  race 

a = 0.8259  mm  , b = 0.0938  mm 
and  for  outer  race, 

a = 0.8722  mm  , b = 0.1131  mm 
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For  Bearing  B at  9.500  rpm 

For  the  operating  conditions: 
axial  preload  : 1,000  N 

bearing  speed  : 9,500  rpm 
operating  contact  angles: 

a.  = 68.73“  , ttg  = 51.70“ 
loads: 

= 247.5  N , Qq  = 300.9  N 
contact  areas: 
for  inner  race 

a = 0.9427  mm  , b = 0.1111  mm 
and  for  outer  race, 

a = 0.9955  mm  , b = 0.1252  mm 

General  Ball  Motion  Based  on  the  Race-Control  Theory 
The  race-control  theory  provides  a convenient  method  to 
compute  the  bearing  kinematics  when  the  gyroscopic  moment  is 
small.  In  this  section,  the  controlling  race  is  assumed  to  be 
the  outer  race.  The  heat  generation  rates  resulting  from  ball 
spinning  are  also  calculated. 

The  rolling,  spinning  (only  on  the  inner  race) , and  the 
spin-to-roll  ratio  for  Bearing  A and  Bearing  B are  shown  in 
Fig. 6-1  and  Fig. 6-2.  Due  to  the  centrifugal  force,  the  spin- 
to-roll  ratios  increase  as  the  speed  increases.  The  heat 
generation  rates  are  shown  in  Fig. 6-3. 
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Figure  6-1  Spinning  and  rolling  velocities  in  Bearing  A (15° 
contact  angle)  and  Bearing  B (60°  contact  angle),  both 
bearings  are  preloaded  by  1000  N. 


Bearing  Speed  RPM  x 1000 


Figure  6-2  Spin-to-roll  ratios  for  the  two  bearings  under 
1000  N preload. 
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Bearing  Speed  RPM  x 1000 

Figure  6-3  Heat  generation  rates  due  to  ball  spinning  of  the 
two  bearings  under  1000  N preload. 
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Time-Simulation  Solution  for  Ball  Skidding  Motion 

As  mentioned,  Bearing  B represents  a bearing  which  failed 
at  9,300  rpm  in  a test  run.  A time  simulation  is  performed  to 
investigate  the  ball  motion  and  heat  generation  rate  in  the 
9,300  rpm  neighborhood. 

Using  Race-Control  Theory  to  Determine  the  Threshold  Speed  of 
Ball  Skidding 

To  determine  the  threshold  speed  of  ball  skidding,  the 
gyroscopic  moment,  J|^  and  the  resistant  torque,  n Q d^, 
in  the  minimum  preload  condition  are  calculated  for  a range 
of  speeds  based  on  the  race-control  theory.  The  minimum 
preload  condition: 

Jfa  0b  “z  ^ M Q db  (2-7) 

For  Bearing  A,  the  results  in  Fig. 6-4  indicate  that  the 
gyroscopic  moment  is  always  smaller  than  the  resistant  torque. 
The  results  for  the  Bearing  B in  Fig.  6-5  indicate  that  the 
gyroscopic  moment  exceeds  the  frictional  torque  at  about  9,000 
rpm,  and  this  means  that  the  ball  skidding  is  expected  to 
occur  at  this  speed. 

Speeds  for  Time  Simulations 

In  order  to  understand  the  ball  motion  at  the  speeds 
around  9,000  rpm,  time  simulations  are  performed  for  the 
following  three  speeds: 
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Figure  6-4  Gyroscopic  moment  and  the  resistant  moment  on 
Bearing  A (15°  contact  angle)  under  1,000  N preload. 


Figure  6-5  Gyroscopic  moment  and  the  resistant  moment  on 
Bearing  B (60°  contact  angle)  under  1,000  N preload. 
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7.500  rpm 

8.500  rpm 

9.500  rpm 

If  there  is  no  radial  load,  the  time-simulation  method 
and  quasi-static  method  are  virtually  the  same  method. 

Format  of  the  Time-Simulation  Output  Plots 

The  time-simulation  outputs  are  plotted  in  a form  that 
can  be  found  in  Fig.  6-6.  The  following  features  of  these  plots 
should  be  noted. 

1. The  horizontal  axis  shows  the  ball  position  which  is  used 

in  place  of  time  (ball  position  is  proportional  to  time) . 
Note  that  the  balls  are  circulating  inside  the  bearing. 

2.  The  initial  conditions,  the  values  at  the  vertical  axis,  are 

data  from  the  race-control  theory. 

3.  The  initial  part  of  each  plot  showing  a transient  variation 

should  be  ignored. 

Time-Simulation  Results  for  no  Radial  Load  Conditions 

Fig. 6-6  shows  that  the  heat  generation  rate  almost 
doubles  (from  32  Watts  to  57  Watts  per  ball)  when  the  speed 
increases  from  8,500  to  9,500  rpm.  This  is  a significant 
increase  of  heat  generation  and  is  believed  to  be  the  cause 
of  the  failure  of  Bearing  B.  It  also  can  be  seen  that  the  heat 
generation  rate  calculated  from  the  time  simulation  is  five 
times  that  calculated  from  the  race-control  theory. 
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The  time  simulation  also  shows  that  the  gyroscopic  moment 
already  increases  the  heat  generation  rate  at  7,500  rpm  which 
is  lower  than  the  speed  predicted  by  the  race-control  theory. 
This  is  because  the  ball  spinning  is  occurring  at  the  same 
time. 

When  the  ball  skidding  occurs,  the  ball  motion 

represented  by  the  four  velocity  components  - Wy,  w„,  and 

X y z 

0^  - is  shown  in  Fig.  6-7.  The  most  important  point  of  this 
figure  is  that  which  represents  the  skidding  speed  is  not 
zero.  Furthermore,  since  other  velocity  components  also  depart 
from  the  solutions  given  by  race-control  theory  (which  implies 
no  gross  sliding) , gross  sliding  is  clearly  seen  to  be  present 
in  this  bearing. 

The  difference  between  the  heat  generation  rates 
calculated  by  the  time-simulation  method  and  the  race-control 
theory  is  shown  in  Fig. 6-8.  This  figure  shows  that  ball 
skidding  causes  substantial  increase  of  heat  generation.  It 
also  can  be  found  that  the  ball  spinning  friction  in  Bearing 
B is  already  higher  than  that  of  Bearing  A. 

Time  Simulation  for  Bearing  B Under  a Radial  Load 

If  a radial  load  is  present,  the  load  distribution  on  the 
Bearing  B is  non-axisymmetrical , and  the  rotational  speeds  of 
each  ball  undergo  a time-cycle  variation. 
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ball  position  (degree  i 1000) 

Figure  6-6  Heat  generation  of  Bearing  B for  different  speeds 
(for  axisymmetrical  loading,  the  time-simulation  and 
quasi-static  method  results  are  the  same) . The 
initial  values  are  the  race-control  theory  results, 
and  the  initial  transient  variation  should  be 
ignored. 


77 


ball  position  (dogree  i 1000) 


Figure  6-7  Time  simulation  solutions  for  ball  rotational 
speeds.  The  shaft  speed  is  9,500  RPM.  The  initial 
values  are  race-control  theory  results,  and  the  initial 
transient  variation  should  be  ignored. 


Heat  generation  rate  per  ba  I I 
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Figure  6-8  Heat  generation  rates  calculated  by  the  race- 
control  theory  and  time  simulation  method.  Both 
bearings  are  under  1000  N axial  preload.  (For  Bearing 
A the  time  simulation  results  are  close  to  the  race- 
control  theory  results) 
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Load  distribution  in  a radially  loaded  bearing.  Fig. 6-9 
and  Fig. 6-10  show  the  loads  and  contact  angles  of  each  ball. 
The  operating  condition  is;  1000  N axial  preload,  500  N radial 
load,  and  9500  rpm.  If  the  radial  load  is  not  present,  the 
loads  and  contact  angles  are 

Q.  = 247.5  N,  Qq  = 300.9  N,  a-  = 68.73°,  = 51.70° 

Heat  generation  in  a radial  loaded  bearing.  The  heat 
generation  in  this  bearing  is  shown  in  Fig. 6-11.  An  enlarged 
section  of  Fig. 6-11  is  shown  in  Fig. 6-12.  It  can  be  seen  that 
the  time  simulation  has  been  performed  for  five  revolutions 
of  the  balls.  The  peaks  of  the  time-simulated  curve  correspond 
to  the  lowest  loaded  position,  and  the  valleys  to  the  highest 
loaded  position. 

Fig. 6-11  shows  that  a radial  load  can  significantly 
increase  the  magnitude  of  the  ball  skidding.  The  maximum  heat 
generation  rate  reaches  100  Watts  (at  the  final  state)  which 
is  twice  the  heat  generated  in  the  pure  axial  load  case. 

The  heat  generation  in  a ball  results  from  three  sources: 
lateral  sliding  and  slidings  in  the  rolling  direction  at  the 
inner  and  outer  races.  This  is  the  reason  why  the  heat 
generation  per  ball  does  not  varies  smoothly  with  the  position 
of  ball  (as  shown  in  Fig. 6-13). 

Comparison  between  the  cniasi-static  solution  and  the 
time-simulation  solution.  In  Fig. 6-11  and  6-12,  the  quasi- 
static solution  is  represented  by  the  curve  marked  with  the 
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"+”  symbol  (each  symbol  represents  a ball) . It  can  be  easily 
seen  that  the  quasi-static  results  are  much  smaller  than  the 
time-simulation  results.  Fig.6-14a  and  b provide  an 
explanation.  The  quasi-static  solution  is  obtained  by 
calculating  the  equilibrium  conditions  for  each  ball  at 
different  positions  and  does  not  consider  the  "history”  of 
ball  motion.  However,  sometimes  (as  in  this  case)  the  friction 
torque  between  balls  and  races  is  not  enough  to  change  the 
ball  spinning  axis,  and  in  such  a situation,  the  variation  of 
the  direction  of  the  ball  spinning  axis  deviates  from  the 
quasi-static  solution. 


Figure  6-9  Load  distribution  on  each  ball  in  Bearing  B (500 
N radial  load,  1000  N preload,  and  9,500  RPM) . If  there 
is  no  radial  load,  = 247.5  N and  = 300.9  N. 
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Figure  6-10  Contact  angles  of  each  ball  in  Bearing  B (500  N 
radial  load,  1000  N preload,  and  9,500  RPM) . If  there 
is  no  radial  load,  a-  = 68.73*  and  a = 51.70*. 


Figure  6-11  Heat  generation  in  Bearing  B when  it  is  under  a 
radial  load  (the  time  simulation  and  quasi-static 
methods  show  different  results) . 
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Figure 


6-12  An  enlarged  section  of  Fig. 6-11. 
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Figure  6-13  Heat  generation  in  a ball  results  from  three 
sources:  lateral  sliding  and  slidings  in  the  rolling 
direction  at  the  inner  race  and  outer  race. 
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a)  quasi -static  solution 


minimum 


450 

3471 

10.32° 
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Figure  6-14  Heat  generation  for  a ball  at  different 
positions  (for  Bearing  B)  : a)  is  calculated  by 
the  quasi-static  method  and  b)  is  calculated 
by  the  time  simulation  method. 
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Temperature  Distribution  in  the  EHD  Film 
The  temperature  distribution  in  the  EHD  contacts  is 
calculated  in  this  section.  The  main  purpose  of  calculating 
temperature  distribution  is  to  determine  the  traction 
coefficient. 


Elastohvdrodvnamic  Film  Thickness 

In  order  to  calculate  the  temperature  distribution,  the 
elastohydrodynamic  film  thickness  has  to  be  computed  first. 
The  following  lubricant  properties  are  used  for  this 
computation: 


T]  (absolute  viscosity  at  atmosphere  condition)  = 0.08 
Pa • sec 

a (viscosity-pressure  factor)  = 2.0  x 10'®  (1/Pa) 


The  Hamrock-Dowson  central  film  thickness  formula  is  used: 

h„ 


= 2.69 


yO.67  qO.53 


R 


W 


0.067 


The  rolling  velocity  can  be  calculated  by  (for  0^^  = 0) 


U = (0ir  - 0b) 


(d_  - D.  cosa«) 


since 


0b  = 


(dp  - Dfa  cosap) 
2 d_ 


^ii 


u = 


V " cosaQ)2  ^ 
0 


4 d. 


ir 
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Due  to  the  high  rolling  speed,  the  films  have  substantial 
thickness.  For  speeds  from  0 to  15,000  rpm,  the  central  film 
thicknesses  are  shown  in  Fig. 6-15. 


Figure  6-15  Film  thicknesses  for  Bearing  A at  different 
speeds  and  under  1000  N preload. 
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Temperature  Distribution  in  the  Lubricating  Film 

The  following  data  are  used  in  the  temperature 
calculation: 

Material  properties  of  ball,  race  and  lubricant: 
for  steel  ball  and  race, 
p = 7800  kg/m^ 

C = 470  J/kg  *C 
kj  = 43  W/m  °C 
for  lubricant, 

p =800  kg/m^ 

C = 2000  J/kg 
k^  = 0.14  W/m 
Coefficient  of  friction: 

The  coefficient  of  friction  is  assumed  to  be  0.054. 

Temperature  Rise  Inside  the  Contact  Caused  bv  Ball  Spinning 
The  following  condition  occurs  in  the  inner  race  contact 
of  Bearing  A when  it  operates  at  10,000  rpm. 

Q = 163.7 

u^  (rolling  speed)  = 29.8  m/sec 
h = 2.6  /im 

Uj  (spinning  speed)  = 938.1  rad/sec 
a = 0.8259  mm,  b = 0.0938  mm. 
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The  ball  spinning  speed  creates  the  sliding  speeds  in  the 
contact  as  shown  in  Fig. 6-16. 

Using  the  above  spinning  speed,  Eg. (4-4) , and  Eg. (4-6) , 
the  temperature  distribution  inside  the  inner  race  EHD  contact 
is  plotted  in  Fig. 6-17  and  Fig. 6-18.  The  maximum  temperature 
occurs  at  X = ± 0.65  a and  y = 0.25  b.  The  temperature  is 
calculated  at  the  mid-plane  of  the  film.  It  also  can  be  found 
that  the  calculated  film  temperature  is  much  higher  than  the 
surface  temperature.  This  is  because  the  conductivity  of  the 
lubricant  is  only  0.3%  of  the  metal. 


(x/a)  along  y = 0 

Figure  6-16  Sliding  speeds  in  a contact  caused  by  ball 
spinning:  u = 938.1  rad/sec,  a = 0.8259  mm  (for 
Bearing  A at  10,000  RPM) . 
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Figure  6-17  Film  temperature  distribution  along 

y = 0.25  b (transverse  direction),  resulting  from  the 
spinning  in  Fig. 6-16  (Bearing  A under  1,000  N preload) . 


Figure  6-18  Film  temperature  distribution  along  x=o.65  a 
(rolling  direction) , resulting  from  the  spinning  in 
Fig. 6-16  (Bearing  A under  1,000  N preload). 
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Traction  Curves  Considering  the  Thermal  Effects 

In  order  to  determine  the  effects  of  sliding  inside  the 
EHD  contact,  two  traction  curves  are  calculated  by  using  the 
following  parameters: 

1)  limiting  shear  stress  = 0.054  P. 

2)  no  limiting  shear  stress 
Lubricant  viscosity  model: 

ri  (viscosity  at  40°  C)  = 0.08  Pa«s 
a (viscosity-pressure  factor)  = 2.0  x 10‘®  (1/Pa) 

)3  (viscosity-temperature  factor)  = 0.04  (1/°C) 

Other  parameters: 
h = 2.6  /im 

U (rolling  speed)  =30  m/sec 
b = 0 . 1 mm 

The  values  of  /)  C k are  the  same  as  that  specified 
earlier. 

The  procedure  used  in  the  calculation  of  traction 
coefficient  is  outlined  in  the  flow  chart  in  Fig. 3-3.  The 
maximum  values  of  the  two  traction  curves  in  Fig. 6-19  occur 
when  the  sliding  speed  is  about  0.2  m/sec.  It  can  be  found 
that  the  limiting  shear  stress  significantly  affects  the 
maximum  value  of  the  traction  coefficient.  The  two  traction 
curves  exhibit  negative  slopes  at  high  sliding  speeds. 
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Effects  of  the  Shape  of  the  Traction  Curve 

To  investigate  how  the  variation  of  the  traction  curve 
affects  the  ball  skidding  phenomenon,  time  simulations  are 
performed  on  the  following  three  assumed  traction  curves  (in 
Fig. 6-20) : 

Curve  1:  = 0.2  m/sec,  /Xq  = 0,  = 0.054,  = 0.054 

Curve  2:  u^^^  = 0.2  m/sec,  /Xq  = 0'  " 0.075, 

= 0.035  (when  sliding  > 3 m/sec) 

Curve  3:  U|^  = 0.2  m/sec,  /Xq  = 0,  /x^  = 0.075, 

/Xf  = 0.035  (when  sliding  > 2 m/sec) 

The  heat  generation  rates  on  Bearing  B at  9,500  rpm 
corresponding  to  the  above  three  traction  curves  are  presented 
in  Fig. 6-21.  It  has  been  shown  in  Fig. 6-6  that  Traction  Curve 
1 will  lead  to  ball  skidding.  An  interesting  result  is  that 
the  minor  difference  between  Traction  Curves  2 and  3 creates 
a large  difference  in  the  heat  generation  rate  (in  Fig. 6-21) . 
This  is  because  only  the  traction  coefficients  in  the  range 
between  0.5  m/sec  to  2 m/sec  are  really  important  to  the  ball 
skidding,  due  to  the  ball  spinning. 

It  also  can  be  found  that  the  negative  slope  in  the 
traction  curve  does  not  have  significant  effects  on  the  ball 
skidding  phenomenon. 
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Figure  6-19  Calculated  traction  curves  (using  Roelands 
viscosity) : in  one  curve,  the  limiting  shear  stress  is 
assumed  to  be  0.54  p;  in  another  curve,  there  is  no 
limiting  shear  stress. 


Figure  6—20  Assumed  three  traction  curves  for  investigating 
how  the  shape  of  traction  curve  affects  ball  skidding. 
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ball  position  (degree  i 1000) 

Figure  6-21  Heat  generation  rates  calculated  according  to  the 
three  traction  curves  in  Fig. 6— 20.  The  minor  difference 
between  Traction  Curves  2 and  3 results  in  a significant 
difference  in  heat  generation  rate.  This  result  shows 
that  the  traction  coefficients  in  the  thermal  region  is 
important  to  the  ball  skidding.  The  initial  transient 
variation  should  be  ignored. 
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Film  Temperature  and  Estimated  Bulk  Temperature 

The  heat  generation  rates  in  Bearing  A and  Bearing  B at 
different  speeds  have  been  shown  in  Fig.  6-8.  It  is  also 
desirable  to  know  the  temperatures  in  EHD  films  resulting  from 
the  calculated  heat  generation  rates. 

Bulk  bearing  temperature  . It  is  extremely  difficult  to 
calculate  the  bulk  bearing  temperature,  thus  this  temperature 
is  estimated  by  assuming  that  all  the  heat  is  taken  away  by 
a lubricant  flow  of  0.1  liter/min.  The  bulk  bearing 
temperature  is  calculated  by 


H 

■^bulk  = 

^oil  *^oil  ® 


+ T 


oil 


(4-1) 


This  estimation  method  gives  bulk  bearing  temperatures 
that  are  close  to  the  true  measured  values  from  bearing  outer 
races. 

Film  temperature.  For  Bearing  A,  the  calculated 
temperature  rise  in  the  lubricating  film  is  40 ”C  at  10,000 
rpm,  and  the  estimated  bulk  temperature  is  30“C  (Fig. 6-22). 
There  shouldn't  be  significant  lubricant  deterioration 
occurring  due  to  this  temperature. 

The  temperature  rise  in  Bearing  B caused  by  ball  skidding 
is  shown  in  Fig. 6-23.  Due  to  the  high  heat  generation  rate  (57 
Watts  per  ball),  very  high  temperatures  (435°C  bulk 
temperature  and  930 'C  film  temperature)  result.  Note  that  the 
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bulk  temperature  of  this  bearing  at  6,000  rpm  is  20 °C  which 
corresponds  to  the  measured  outer  race  temperature. 


Bearing  Speed  RPM  i 1000 


Figure  6-22  Temperature  rises  in  Bearing  A.  Heat  is  carried 
away  by  an  0.1  lit/min  oil  flow. 


temperature  (lOOG  C) 
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Figure  6-23  Temperature  rises  in  Bearing  B. 
away  by  an  0.1  lit/min  oil  flow. 


Heat  is  carried 


CHAPTER  7 

SUMMARY  AND  CONCLUSIONS 


Based  on  a bearing  dynamic  analysis,  the  ball  motion  and 
heat  generation  rates  in  a 60'  and  a 15"  contact  angle  ball 
bearings  are  calculated.  The  major  purpose  of  this 
investigation  is  to  explain  the  sudden  failure  of  the  60" 
contact  angle  ball  bearing  in  a machine  tool  spindle.  Some 
important  results  are: 

l.The  time-simulation  results  on  Bearing  B (60'  contact  angle) 
show  that  when  the  bearing  is  accelerated  from  the  8,500  rpm 
to  9,500  rpm,  the  heat  generation  rate  almost  doubles  - from 
32  Watts  to  57  Watts  per  ball.  According  to  the  calculated 
heat  generation  rates,  if  the  bulk  temperature  of  Bearing 
B at  6,000  rpm  is  20 "C  (in  fact,  this  is  the  temperature 
measured  at  the  Bearing  B outer  race) , the  bulk  temperature 
at  9,500  rpm  with  the  same  lubricant  flow  rate  will  be 
435 "C.  If  the  bearing  is  non-axisymmetrically  loaded,  the 
temperature  will  be  even  higher. 

It  is  hard  to  pinpoint  at  which  temperature  the  seizure 
will  occur.  Since  the  AISI  M series  steels  (the  common 
bearing  materials)  begin  to  soften  at  400 'C,  it  can  be 
expected  that  large  deformation  of  the  balls  can  be  caused 
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by  the  compressions  between  the  balls  and  races  when  the 
temperature  is  above  this  value. 

2 . It  can  be  found  that  the  ball  spinning  decreases  the 
threshold  speed  of  ball  skidding  in  two  ways.  First,  the 
ball  spinning  takes  away  some  of  the  available  traction 
force  at  the  inner  race  contact.  This  is  the  reason  why  ball 
skidding  actually  occurs  in  Bearing  B at  speeds  lower  than 
9,000  rpm  which  was  predicted  without  considering  the 
combined  effect  of  the  ball  spinning  and  ball  skidding. 
Second,  the  heat  generation  caused  by  the  ball  spinning 
diminishes  the  traction  coefficient.  It  was  shown  that  the 
magnitude  of  the  traction  curve  in  the  thermal  region  is 
important  to  the  threshold  speed  of  ball  skidding.  It  was 
also  found  that  the  negative  slope  of  the  traction  which  may 
sometimes  cause  self-excited  motion  does  not  affect  the  ball 
skidding  phenomenon. 

3.  The  ball  skidding  generates  much  more  heat  than  the  ball 
spinning  does.  The  heat  generation  rates  in  Bearing  B 
calculated  by  time  simulations  are  about  five  times  that 
calculated  by  race-control  theory  which  only  considers  the 
ball  spinning.  The  heat  generation  rate  in  Bearing  B (60° 
ball  bearing)  is  also  much  higher  than  that  in  Bearing  A (57 
Watts  per  ball  at  9,500  rpm  in  Bearing  B compared  to  2.62 
Watts  heat  per  ball  at  10,000  rpm  in  Bearing  A ). 

4.  From  a thermal  analysis,  the  temperature  in  the  EHD  film 
caused  by  the  2.62  Watts  heat  rate  is  40°C  above  the  metal 
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surface  temperature  (rolling  speed  is  30  m/sec) . There  is 
only  a small  increase  (0.8°C)  of  the  temperature  on  the 
metal  surface  in  the  contact  area.  The  difference  is  due  to 
the  fact  that  lubricant  thermal  conductivity  is  usually  only 
0.3%  of  that  of  metal. 

5. Ball  skidding  phenomenon  is  a typical  example  showing  that 
traction  coefficient  is  related  to  the  stability  of  motion. 
In  fact,  there  are  many  instances  in  the  physical  world 
(eg. , skidding  at  the  traction-driving  contacts,  stick-slip 
motion  [30],  self-excited  cage  whirl  motion  [13,14],  etc.) 
that  display  this  characteristic.  Thus,  the  mathematical 
modelling  technique  used  in  this  work  should  find 
applications  in  areas  other  than  ball  bearings. 


APPENDIX  A 

HERTZIAN  POINT-CONTACT  SOLUTIONS 


This  appendix  presents  the  Hertzian  solutions  for  the 
contact  geometry  and  pressure  distribution  between  two 
contacting  elliptic  bodies. 


b)  equivalent  contact 


Figure  A-1  Contact  geometry  for:  a)  an  original  contact;  b) 
the  equivalent  contact. 
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Equivalent  Radius  for  Two  Elliptic  Bodies  in  Contact 

When  two  elliptical  bodies  are  brought  together  into 
contact  as  shown  in  Fig.A-la,  the  contact  geometry  can  be 
represented  by  an  equivalent  elliptical  body  near  a plane  as 
in  Fig.A-lb.  The  geometric  separation  between  the  two  original 
elliptical  bodies  is  the  same  as  the  separation  between  the 
equivalent  body  and  the  plane.  At  a point  (x,y)  , the  geometric 
separation  can  be  written  as 


2 2 

X y^ 


S(x,y)  - + 

(A-1) 

2 R, 

2 Ry 

In  deriving  the 

above  expression,  the 

following 

approximation  has  been  employed: 

R,  - y R,2  - . 

x2 

Rx  ( Rj(  + . . 

2 R^ 

• ) 

Since  the  contact 

region  is  usually  small, 

the  above 

approximation  has  sufficient  accuracy.  A similar  approximation 

for  Ry  has  also  been  used. 

The  radius  of  the 

equivalent  elliptical  body  can  be 

calculated  by 

1 1 

1 

(A-2a) 

+ 

^xb 

r* 

xr 

1 1 

1 

(A-2b) 

cr 

r 

yr 

1 1 

1 

(A-2c) 

+ 

R R^ 

R, 
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where  R^,  are  the  equivalent  radii  in  the  x-direction  and 
y-direction  respectively 
R is  the  over  all  equivalent  radius 

^bx  ^by  radii  of  the  ball  in  the  planes 

containing  x-axis  and  y-axis  respectively, 
r^^  and  r^y  are  the  radii  of  race  curvature  in  the  planes 
containing  x-axis  and  y-axis  respectively. 


Figure  A-2  Deformation  on  a sphere  when  it  is  pressed  into  a 
semi-infinite  surface. 

Elasticity  Theory 

For  calculating  the  elastic  deformation  caused  by  an 
applied  load,  the  following  equation  presented  by  Timoshenko 
and  Goodier  [29]  can  be  applied.  This  equation  specifies  the 
deformation  at  a point  on  a semi-infinite  surface  caused  by 
a pressure  p acting  on  a small  area  dA  from  a distance  r. 


2 p (1-1/2) 
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dS 


or 


7T  E'r 


dA 


(A-3) 


2 p dxi  dyj 

J J , i 5 (A-4) 

IT  E'  J (y-yj)'^  + (X-Xj)^ 


where  E'  is  defined  as  E/(l-u^),  the  elastic  modulus,  and  u is 
the  Poison's  ratio. 


Hertzian  Pressure  Distribution 

As  shown  in  Fig.A-2,  an  elliptical  body  is  pressed  into 
a semi-infinite  surface  by  a distance  S.  The  deformation  of 
a point  at  a distance  r from  the  center  point  can  be  expressed 
as 

2 2 
X y 

Wj  + W2  = 6 (A-5) 

2 2 Ry 

where  Wj  and  W2  represent  the  deformation  on  the  elliptical 

body  and  the  plane  respectively.  It  was  found  by  using 

Eq. (A-5)  that  the  pressure  distribution  which  engenders  such 

deformation  is 

3Q  X y 

a = [ 1 - ( )2  - ( )2  jl/2 

27rab  a b 


Then,  the  maximum  pressure  of  the  contact  is 
3Q 

p 

max 


2jrab 


(A-7) 


103 


where  a is  the  semi-major  axis  of  the  contact  ellipse,  and  b 
is  the  semi-minor  axis  of  the  contact  ellipse. 

Hertzian  Solution  for  the  Point-Contact  Application 

Harris  [l]  has  summarized  the  solutions  for  the  contact 
and  deflection  for  the  application  of  ball  bearings  as 
follows: 


a = ( ) a*  ( Q R ) 

E' 

(A-8) 

3 

b = ( ) b*  ( Q R ) 

E' 

(A-9) 

elastic  deformation 

3 , q2/3 

S - ( s* 

E'  2 r1/3 

(A-10) 

Eq. (A-10)  can  be  rewritten  in  the 

form  of  Eq. (3-2)  as 

Q = K 5^/^ 

E'  2 R^/^ 

and  K - ( ) 

3 6*3/2 

(A-11) 

where  E'  is  defined  as: 
E 


E'  = 


1 - u" 


(A-12) 


The  constants  a*,  b*,  and  S*  are  defined  as 
a*  = ( 


2 E 


2 jl/3 


(A-13a) 


TT 
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b* 


S* 


2 E, 


1/3 


7T  K 
2 Ej 
7T 


2 El 


1/3 


(A-13b) 

(A-13C) 


Where  R is  the  equivalent  radius  of  the  elliptic  contact  as 
defined  earlier,  and  k is  the  ellipticity  ratio  defined  as 


a 

K = (A-14) 

b 


Ej  is  the  complete  elliptical  integration  of  the  first 

kind 


1 

Ej  = J —7^--  5 d0  (A-15a) 

y [1-  ( 1-  bva^  ) sin‘^0] 


and  E2  is  the  complete  elliptical  integration  of  the  second 
kind 

E2  = J y [1-  ( 1 -h^/a^  ) sin^0]  d0  (A-15b) 


In  the  above  equations,  the  ratio  k is  an  unknown  which 
has  to  be  determined  by  solving  the  following  equation: 

( li^  + 1 ) E2  - 2 Ej 

y = F(R)  (A-16) 

( K - 1 ) E2 


where,  F(R)  is  defined  as 

1 1 

^bx  ^by 


1 


1 


rx 


ry 


F(R)  = [( 


) + ( 


)]  R 
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Curve-Fitting  Formulas  for  Elliptic-Contact  Deformation: 

Harris  [1]  has  solved  Eq. (A-16) , and  he  presented  the 
solutions  in  tables.  A simpler  method  to  obtain  the  values  of 
K,  Ej,  and  E2  is  to  use  the  equations  derived  by  Brewe  and 
Hamrock  (1977)  as  shown  below.  These  equations  are  obtained 
by  using  curve-fitting  with  the  least  square  method. 


R, 


X ,0.6360 


K,  = 1.0339  ( 


(A-17a) 


R, 


‘y 


E,  = 1.5277  + 0.6023  ln(Ry  / R^) 
•i  X y 


(A-17b) 


0.5968 


E2  = 1.0003  + 


(A-17C) 


(V^y) 


APPENDIX  B 

EHD  LUBRICATION  THEORY 


Since  the  mid-20  century,  the  elastohydrodynamic  lubrica- 
tion phenomenon  has  been  intensively  studied.  It  is  now 
commonly  known  that  the  development  of  the  EHD  film  is 
governed  by  the  following  equations; 

.Lubrication  action  (governed  by  the  Reynolds  equation) 
.Elastic  theory 

.Lubricant  Rheology  (variation  of  lubricant  viscosity) 
.Energy  equation  (determines  the  temperature  distribution 
and  the  viscosity  in  the  film) 

The  Reynolds  equation  is  derived  by  combining  the  Navier- 
Stokes  equation  and  the  continuity  equation.  In  deriving  the 
Reynolds  equation,  the  fluid  is  assumed  to  be  Newtonian  and 
iso-viscous.  The  Navier-Stokes  equations  for  a Newtonian  flow 
are  written  as 

3p  d 3u 


= ( T]  ) 

5x  dz  dz 

dp  d dv 

= ( T]  ) 

dy  dz  dz 


The  Reynolds  equation  is  derived  by  using  the  integral 
form  of  the  continuity  equation; 


106 


107 


d d d 

f u dz  + f V dz  + r w dz  = 0 

dx  dy  dz 

Also,  according  to  Fig.B-1,  the  boundary  conditions  are 

z = 0,  u = Uj,  V = Vj,  w = Wj 
z=h,  u = U2,  v = V2,  w = W2 


ball 


Figure  B-1  Coordinate  System  for  deriving  the  Reynolds 
equation. 

Derivation  of  the  Reynolds  Ecmation 

According  to  the  Navier-Stokes  equation,  the  velocity 

components  are  in  parabolic  forms  as 

1 3p  1 

— + Cj  z + <^2 


u 


2 r?  dx 
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1 3p 

2 r?  ay 


+ c- 


1 


+ c 
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By  utilizing  the  boundary  conditions,  the  constants  can 
be  determined.  The  velocity  profiles  now  can  be  written  as 


h - z 3p 

u = - z ( ) 

2 T]  dx 

h - z ap 

V = - z ( ) 

2 T]  ay 


h - z z 

+ Uj  ( ) + U2  


h - z 

+ Vj  ( ) + V2 


h 

z 


By  substituting  the  above  velocity  profiles  into  the 
integral  form  of  continuity  equation,  we  arrive  at  the 
standard  form  of  Reynolds  equation: 


3 ap  a ap  ah 

( h-^  ) + ( h^  ) = 12  (Uj  + U2)  + 

dx  T}  dx  dy  r]  dy  dx 

dh 

12  (Vj  + V2)  + (Wj  - W2)  (B-1) 

ay 


Solution  Methods  for  the  Reynolds  Equation 

Solving  the  Reynolds  equation  for  the  EHD  lubrication  is 
difficult  due  to  the  following  two  reasons:  1.  a small  change 
in  the  film  shape  creates  a significant  change  in  the  pressure 
distribution.  2.  lubricant  viscosity  increases  with  pressure. 

Thus,  the  following  procedure  which  is  effective  for  HD 
problems  but  is  not  applicable  in  EHD  problems. 
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film  shape  ->  pressure  ->  deformation  -> 
film  shape  ->  ... 

The  above  method  is  known  as  the  forward  method  which 
usually  does  not  lead  to  a convergent  solution.  It  has  been 
suggested  that  an  inverse  method  should  be  used  in  order  to 
avoid  the  divergence  problem.  The  inverse  method  utilizes  the 
Reynolds  equation  to  determine  film  shapes  from  pressure 
profiles  instead  of  obtaining  pressure  profiles  from  film 
shapes  as  in  the  forward  method. 

Film  Shape 

The  shape  of  an  EHD  film  can  be  written  as 

h(x,y)  = hQ  + 5(x,y)  + S(x,y)  (B-2) 

where,  hg  is  a constant  to  be  determined. 

S(x,y)  is  the  elastic  deformation 

S(x,y)  is  the  initial  separation  which  can  be  written 
as 

2 2 
X y 

S(x,y)  = + 

2R  2R 


Determination  of  Film  Shape  bv  the  Inverse  Method 

The  purpose  of  using  the  inverse  method  is  to  obtain  a 
film  shape  from  a prescribed  pressure  profile.  The  inverse 
method  needs  to  use  the  integrated  Reynolds  equation  (one- 
dimensional) as, 
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ap 


ay 


12  rj  V 


(B-3) 


where  hjjj  is  the  film  thickness  at  dp/dy  = 0 (maximum  pressure 
point) . The  integrated  Reynolds  equation  is  derived  by 
performing  an  integration  over  y. 

From  the  film  shape  equation,  Eq. (B-2) , it  can  be  found 
that  when  the  pressure  center  (the  deformation  center)  moves 
forward,  the  slope  of  the  film  shape  changes.  This  fact  can 
be  utilized  in  the  selection  of  the  pressure  profile  in  the 
analysis  of  the  inverse  method. 
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